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Abstract: Belinskii, Khalatnikov, and Lifshitz (BKL) conjectured 
that the description of the asymptotic behavior of a generic solution 
of Einstein equations near a spacehke singularity could be drastically 
simplified by considering that the time derivatives of the metric asymp- 
totically dominate (except at a sequence of instants, in the 'chaotic case') 
over the spatial derivatives. We present a precise formulation of the BKL 
conjecture (in the chaotic case) that consists of basically three elements: 
(i) we parametrize the spatial metric gij by means of Iwasawa variables 
(j3"',M"'i); (ii) we define, at each spatial point, a (chaotic) asymptotic 
evolution system made of ordinary differential equations for the Iwa- 
sawa variables; and (iii) we characterize the exact Einstein solutions 
P, J\f whose asymptotic behavior is described by a solution /3[o], A/jo] of 
the previous evolution system by means of a ^generalized Fuchsian sys- 
tem^ for the differenced variables (3 = (3 — /?[o] , M = M — M[o] , and by re- 
quiring that P and M tend to zero on the singularity. We also show that, 
in spite of the apparently chaotic infinite succession of 'Kasner epochs' 
near the singularity, there exists a well-defined asymptotic geometrical 
structure on the singularity : it is described by a partially framed flag. 
Our treatment encompasses Einstein-matter systems (comprising scalar 
and p-forms), and also shows how the use of Iwasawa variables can sim- 
plify the usual ('asymptotically velocity term dominated') description of 
non-chaotic systems. 
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1 Introduction 



The works [H El E] of Belinskii, Khalatnikov and Lifshitz (BKL) proposed a description of 
the asymptotic behavior of the gravitational field in the vicinity of a spacelike singularity 
of a D = 4 spacetime satisfying the vacuum Einstein equations. They also investigated the 
D = 5 vacuum Einstein case as well as the D = 4 spacetime with a massless scalar field [3] . 
Finally, they analyzed more general Einstein-matter systems (e.g. Einstein- Yang-Mills) in 
[5] . Extension of the BKL analysis to higher dimensions was addressed within the context 
of pure gravity in [6l [7] . It is convenient to express the BKL analysis using Hamiltonian 
methods. This leads to considering the motion of a particle in an auxiliary Lorentzian 
space submitted to the influence of a linear superposition of exponential potential walls 
[B El El Uni EH H] ■ This approach allows for a relatively easy generalization of the work of 
BKL to any spacetime dimension and with any p-form field content [12] . 

As argued by BKL, a drastic simplification in the Einstein equations occurs near a spacelike 
singularity (located at proper time t = 0) in that the partial differential equations for the 
metric can be essentially replaced by ordinary differential equations with respect to time. In 
physical terms this corresponds to an effective decoupling of spatial points xi ^ X2 as t ^ 0. 
Depending on the specific theory at hand (spacetime dimension, field content, couplings to 
the dilatons), the BKL approach leads to expect two possible types of behavior: 

• 'Non-chaotic behavior' (or Monotonic power law): the spatial scale factors (and the 
dilaton fields e'^ if any), behave at each spatial point in a monotone, power-law fash- 
ion in terms of the proper time as one approaches the singularity at t = 0, i.e. at 
each spatial point the asymptotic form of the metric looks like a Kasner metric. On 
the other hand, the p-form fields A have limits as t ^ 0. Theories exhibiting this 
asymptotic behavior are, for instance, pure gravity in D > 11 [6l[71[T3], and gravity 
coupled to a scalar field in any dimensions [H [H] . 

• 'Chaotic behavior': at each spatial point, the asymptotic behavior is given by a chaotic 
succession of an infinite number of increasingly shorter Kasner regimes as one goes to 
the singularity. Important examples of theories exhibiting this asymptotic behavior 
are pure gravity in < 10 [6l[7], and the bosonic sector of all super gravities associated 
with the low energy limit of string or M-theory |15j . 

The non-chaotic case has been formulated in rigorous mathematical terms by consider- 
ing an auxiliary asymptotic dynamics called the 'asymptotically velocity term dominated' 
(AVTD) system |16) . The AVTD system is obtained by neglecting all the spatial deriva- 
tives in the considered Einstein-matter system. Einstein equations then reduce to ordinary 
differential equations (ODEs). The solutions of this asymptotic system are precisely given 
by Kasner-like metrics. Fuchsian methods [T7] can then be used to prove that, given a 
solution of the velocity dominated system, there exists a (geometrically unique) solution of 
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Einstein's equations that asymptotically approaches this solution. These Fuchsian methods 
have been used to mathematically describe cosmological singularities in various simplified 
contexts: Gowdy spacetimes |17] . plane symmetric spacetimes with a massless scalar field 
[H], polarized and half-polarized U{1) symmetric vacuum spacetimes [19], spacetimes with 
collisionless matter and spherical, plane or hyperbolic symmetry [20) . and a particular sub- 
set of general Gowdy spacetimes [21]. It has also been possible to use Fuchsian methods 
to mathematically describe singularities without any symmetries: notably for the Einstein- 
scalar system [T3j, and for many Einstein-matter models including pure gravity in D > 11 
dimensions [13] . 

By contrast, the general inhomogeneous chaotic case has not yet been tackled by rigorous 
mathematical methods. The BKL conjectural behavior has been consolidated by recent 
Iwasawa-variable based analytical treatments [12] and is also supported by numerous nu- 
merical results |22 p23ll2^ l25 t [2 ^ l27|. However there exist neither a clear general formulation 
of the precise asymptotic behavior advocated in the BKL approach, nor any mathematical 
theorems concerning its compatibility with Einstein field equations. The purpose of this pa- 
per is to present a mathematically precise formulation of the BKL conjecture in the chaotic 
case. In other words, we aim at providing a chaotic analog of the AVTD formulation of the 
non-chaotic case. More precisely, we shall describe the asymptotic dynamics of the gravita- 
tional field for Einstein-matter systems, at each spatial point, by a well defined asymptotic 
evolution system made of ODEs . 

In addition to formulating a precise conjecture for the chaotic BKL behavior, we also ad- 
dress the question of whether or not geometrical structures can be defined at the singularity 
and what are these asymptotic geometrical structures. According to the billiard picture, 
most of the metric variables possess well defined limits at the singularity: the ^ off- diagonal 
variables', i.e. all the variables except the diagonal metric components (and the dilatonj^. 
This means in particular that for these variables, initial data can be assigned at the sin- 
gularity. The other variables, i.e. the diagonal variables, have no limit at the singularity. 
Although the 'off-diagonal' variables have finite limits, they are (co)frame (and gauge) de- 
pendent and thereby they do not have a priori a clear geometrical meaning. Nevertheless, 
we can wonder whether it is possible to extract some geometrical information from these 
asymptotic values. It turns out that this is possible, but that this asymptotic geometrical 
structure is less 'rigid' in the chaotic case at hand than it was in the non-chaotic case. 
In the non-chaotic case, the asymptotic geometrical structure is simple to describe. The 
solution is asymptotically given, at each spatial point, by a Kasner-like metric |14tll3j. The 
(spatial) Kasner metric is, in d spatial dimensions, 

gij{t) = t^PHilj + t^P^miTUj + ... + t'^P'^rirj , (1) 
short review of the biUiard picture is presented in section [2] 
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where the pj's {i = 1, 2, 3) are the Kasner exponents subject to the Kasner conditions. [Note 
that this metric possesses a curvature singularity at t = and that the distances are no 
longer defined at this singularity, since either gij ^ cxd or gij 0.] The Kasner coframes, 
i.e. the coframes that diagonalise, at each spatial point, the second fundamental form kij 
with respect to gij have finite limits at the singularity (up to independent rescalings they 
are simply given by u)]^ = Udx"^ , ioj^ = midx"^, = ridx'^) and therefore provide a basis 

of preferred directions, i.e. a 'directional frame' (and co-frame). See Figure [TJ 



(a) (b) 




Figure 1: (a) Non-chaotic behavior : sufficiently close to the singularity, the dynamics of the 
gravitational field can be approximated by a Kasner-like metric at each spatial point. Let us focus 
on one particular spatial point where asymptotically the metric is given by ds^p^^^^i = {t^'P^Ulj + 
t^P^minij + ... + t^'P'^rirj)dx^ dxK When i ^ 0, the directions for which the Kasner exponent pi 
is negative are stretched while the ones with positive exponent are squeezed. At the singularity, 
these directions are still defined, (b) Chaotic behavior : now instead of a Kasner-like metric at 
each spatial point, there is a never ending chaotic succession of Kasner epochs before reaching the 
singularity. Are there still some preferred directions at the singularity? or is some other structure 
of the metric preserved asymptotically? 

For chaotic systems, we could have expected, from the BKL description of the asymptotic 
dynamics of the metric as a never ending chaotic succession of Kasner epochs at each spatial 
point, that no privileged directions can be defined at the singularity especially in view of 
the effect, discovered in p], of a 'rotation' of Kasner frames between two successive Kasner 
epochs. However, we show in section O that an asymptotic geometrical structure can be 
defined at the singularity. This structure is less precise than a frame but more precise than 
a flag and therefore we will call it a partially framed flag. The precise meaning of this notion 
is explained in the sequel. 

This paper is organised as follows. We first review in section [2] the Iwasawa-variables 'cos- 
mological billiards' of [12] in order to introduce our notation and stress important features 
for our purposes. Then, to gain some intuition for how to define the asymptotic system 
of evolution equations parametrizing a 'generic' solution of an Einstein-matter system in 
the chaotic case, we revisit in section [3] the non-chaotic case treated in \14\ I13j . Our 
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new approach is based on Hamiltonian methods and Iwasawa variables, which simplify the 



equations and we rewrite the Hamiltonian Einstein-matter evolution equations in terms of 
the difference between the solution of the full evolution equations and the solution of this 
asymptotic system. Then we argue that the so-obtained 'differenced system' of equations 
is of the Fuchsian typ^. We also treat the constraints by defining asymptotic constraints 
which, when they are satisfied, imply the vanishing of the exact constraints. Next we turn 
to our main purpose in section [H that is to give a mathematically precise formulation of 
the chaotic BKL behavior. This is achieved by defining again an asymptotic system of 
evolution equations which is a system of ordinary differential equations though it is not 
necessarily an AVTD system. We then formally rewrite the Hamilton equations in terms of 
the difference between the solution of the full evolution equations and the solution of the 
asymptotic system. Finally, we argue that a stronger version of the usual Fuchs theorem 
is likely to remain valid in the chaotic case. In the last section, we show that for chaotic 
systems, partially framed flags are the asymptotic geometrical structures that stay well de- 
fined at the singularity. 

We mention the recent paper [25] that establishes the relationship between the Iwasawa- 
based billiard approach used here and the dynamical systems approach to cosmological 
singularities. 

2 Appearance of Toda— like walls in Einstein— matter Hamil- 
tonians in Iwasawa variables 

The general systems considered are of the following form 



Units are chosen such that IGttGn = 1, Gn is Newton's constant and the spacetime dimen- 
sion D = d + 1 IS left unspecified. Besides the standard Einstein-Hilbert term the above 
Lagrangian contains a dilaton field (p and a number of p-form fields B'j^l...^^ (for p > 0). 
The generalization to any number of dilatons is straightforward. The p-form field strengths 
iJ'fp) = dB^'^^ are normalised as 



previously 




precisely, we define an asymptotic system of evolution 






pip) 



{p+l)d[^^B 




lb p permutations . 



comparison between the two analyses is done in appendix [Cl 
^We recall Fuchs' theorem in appendix [B] 
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As a convenient common formulation we have adopted the Einstein conformal frame and 
normahsed the kinetic term of the dilaton (p with weight one with respect to the Ricci scalar. 
The Einstein metric g^^i, has Lorentz signature ( — !-•••+) and is used to lower or raise the 
indices; its determinant is denoted by The dots in the action ([2]) above indicate 

possible modifications of the field strength by additional Yang-Mills or Chapline-Manton- 
type couplings [23 [3D] . The real parameter Xp measures the strength of the coupling of 
B'-''^ to the dilaton. In the following, for simplicity, we shall treat the case where there is 
no dilaton (p and indicate what changes occur when (/> is present. 

2.1 Iwasawa variables 

Let us give a schematic review of the Iwasawa-variable cosmological billiards. For a detailed 
derivation, we refer to [12]. 

We choose a slicing of the spacetime we want to construct, A4d = x M, such that the 
singularity occurs at the coordinate time r = +oo. We shall define the time slicing r by 
requiring that the 'rescaled lapse' N = N/^ (where g =detgij) is equal to some given 
(weight -1) time-independent density ;U_i(x) on M^. For simplicity, we take fj,-i{x) = 1 
in the coframe we use, so that = 1 <^=^> N = In other words our time 

coordinate is linked to the 'proper time' dt = —Ndr by dr = —dt/y/g. The slicing is 
built by use of pseudo-Gaussian coordinates defined by a vanishing shift A* = 0, lapse 
A(r, X*) = -y/ g{T, a;*)^_i(x*) = \/ g{T, x*) and metric 

ds^ = -{N{T,x')dTf + g,j{T,x')oj\x^)J{x^) . (3) 

Here w*(x) = uj''j{x)dx^ is a coframe on the given (analytic) spatial manifold M^Q One of 
the useful technical tools we shall employ here consists in replacing the d{d + l)/2 metric 
variables gij, by a new set of variables: d 'diagonal degrees of freedom' P"", together with 
d{d — l)/2 'off-diagonal degrees of freedom' M°'i where M is restricted to be an upper 
triangular matrix (AA'^j = 0, if i < a) with ones on the diagonal {M"'i = 1, if a = i), such 
that 

5.. = Ee-^^"AA^AA'^,. (4) 

a=l 

We shall refer to the algebraic decomposition (jj]) as the Iwasawa decomposition of the 

*Note that in [T2], a coordinate basis is used instead of a general basis lj* but the generahzation is 
straightforward as long as tJ* = lo' j{x'')dx^ does not depend on time. 
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metricij. In d = 3, the components of the metric read exphcitly, 
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^2^-2/3l 

from which one gets (uniquely) 



■2^^|2^-2/32 ^ 



(5) 
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. r2 523511 - 512513 

A/ 3 = 2 — • 

511522 - 5l2 



(6) 



In the cosmological context, one could refer to the 'diagonal metric variables' e~^" as the 
'scale factors'. In |12j the /3's and the dilaton (when present) were collectively denoted 
= {P^, (/>}. In the case considered here (no dilaton), we shall use labels from the 
beginning of the latin alphabet (a, b, c, .., e) to denote the 'diagonal variables' P'^. All other 
variables are called 'off-diagonal variables' and are denoted Q, 

Q = {M, B'"'}, 

where ^S'*"' are the -B'*"' expressed in the generalized Iwasawa coframe 6f^^, 9^^^ := M'^jU)^, 
e.g. we have Bi^, i^ =: M"'^i^...M"'^ i^Bai,,,ap- Note that, in the Iwasawa coframe, the metric 
is diagonal: g""^,^ = e'"^!^" 5ab- 



2.2 Hamiltonian approach in Iwasawa variables 

The Hamiltonian action corresponding to the action ([2]) in any pseudo-Gaussian gauge, and 
in the temporal gauge for the form fields (-Bon...j„_i = 0), readqj 



1 



ji---jp t>(p) 

h---'jp 



I (p) 



H 



(7) 



^Indeed, it is linked to the Iwasawa decomposition of the vielbeins V"i £ SL{d,W) (such that gij = 
ViV"j) which reads V = KAN where K £ SO{d), is a diagonal matrix and A*' is a nilpotent matrix. The 
Iwasawa variables are uniquely specified by requiring that K £ SO(d) be the unit matrix. 

®The term should be added (inside the parenthesis) in the action if a dilaton is present. 
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where the Hamiltonian densitj0 H 



H = NH (8) 

p 

- gR + y g ,^ Ff\ F^^) J'^-^+i (9) 
^ 2 (p+ 1)! Ji-Jp+i 

where i? is the spatial curvature scalar. The dynamical equations of motion are obtained 
by varying the above action w.r.t. the spatial metric components, the spatial p-form com- 
ponents and their conjugate momenta. In addition, there are constraints on the dynamical 
variables, 

~ ("Hamiltonian constraint"), (10) 
"Hj ("momentum constraint"), (11) 
'^(p) ^""' w ("Gauss law" for each p-form) (12) 



wit 



:= -2.^,,+E^<)-^^4^... (13) 



P 

p 



where the subscript \j stands for the spatially covariant derivative. 



As shown in [12] (and as we will see explicitly for some of the terms below) the Hamiltonian 
density of weight 2, 7i, expressed in the Iwasawa variables has the following structure: 

n[P,Q;TT,P]=IC + V, (15) 

where 

V = j;cA(g,P,5./?,92/3,9Q,52Q)e-2--^(^). 

A 

Here G"^ is the inverse of the quadratic form Gab which is defined by Gabd[3"'d(3 := 
Yl'a=i{d'f3"')'^ — (X]a=i '^Z^'^)^- Note the important fact that this metric has a Lorentzian 
signature (—,+,...,+). P stands for {V^a, ^^^^ "^'j, where the P*a's are the momentum 

'^Note that H is & density of weight 1 while 7i is a density of weight 2. 

*If a dilaton is present, the term ^tt^ should be added in K. as well as exponential coupling e'^'''^ in front 
of the term 9 F'^^I.^^^^F'^^'' ^^"'^'•+^ in M and e'^p"* in front oi i^l'p)'^'' T^i^) j^-j^ in /C.The term gg'^ di(j>dj(t) 
should also be added to M. 

^If there are dilatons, the term -K^dicf) should be added to TLi 
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conjugate to the AA"j's and the £^|^y'"'"s are the tt^^^"*'' — i.e. the momentum conjugate 
to the -Bj-^' ^^'s — expressed in the Iwasawa frame. Note that the V'^s are strictly lower 
diagonal (they exist only when i > a and vanish when i < a). 

Note the special structure of the (weight 2) Hamiltonian density 7i with (i) a kinetic term 
vr^ for a 'point particle' of coordinates (i"" moving in (ii) a sum of 'exponential walls' (for 
their 13 dependence). We shall refer to them as Toda walls e~'^'^^^^^ where the wa{I3) are 
certain linear forms in the /3's. 'Toda' refers to the well known Toda models involving such 
exponential walls. For instance, the kinetic terms of the off-diagonal degrees of freedom 
J\f in H give, when expressed in Iwasawa variables, terms proportional to e"^*^^*""^"-* for 
b > a with coefficient proportional to V^. Therefore the kinetic terms for the AA's furnish 
the walls wsab = — called 'symmetry walls'. The kinetic terms of the p-forms in H 
yield a sum of terms proportional to e^"^"^"^^^ where 'Weai...apiP) = P""^ + ••• + ('electric 
p-form walls'). The curvature term —gR in H gives terms proportional to e~^"'°'"=('^) where 
Wabc{P) = + Y2ej^b c{bj!=c) 1^^ ('curvature walls') and their coefficients, when a ^ h, a ^ c 
and b ^ c, are given by (Ci^a^bc)^- Here, the Ci„a are the structure functions of the Iwasawa 
coframe Bf^^ = M^jO^, d9^^^ = -\ Q^^^bcOt^. ^ ^t.- Note that the structure functions Q^, 
depend on the A^'s and the dxM^s. 



A heuristic analysis of the BKL limit indicates the crucial role played by the linear forms 
WAil3)- Indeed, in this limit the walls become infinitely sharp and are located at the 
hyperplanes given by the linear forms WAiP) = 0, the motion is then restricted to the 
region of /3~space defined by the inequalities {wAiP) > 0}. The set of dominant walls is 
the minimal subset {w^(/?)} — the indices A belong to a subset of the indices A — such 
that the subset of inequalities tf^(/3) > implies the full set of inequalities wa{P) > Vj40 
A crucial consistency condition for these definitions, which is found to be satisfied for all 
models, is that the coefficients of the dominant walls be positive: c_a. > O.We decompose 
the set of indices {A} for the walls into the set of indices for the dominant walls {A} and 
the remaining ones ('subdominant walls') {-4'}. Moreover, another crucial structure of the 
potential is that the dependence of the wall coefficients on spatial derivatives in such that, 

V = j;c^(Q,P,9,.Q)e-2-^(^) + Y,CA'{Q,P,d,P,dl(3,dQ,d^Q)e-^^^'^''^ , (16) 

A A' 

where the coefficients of the dominant walls are found never to depend on {dxP, d'^P}- 

For instance, in the case of pure gravity in d space dimensions, the d dominant walls comprise 

(i) d — 1 dominant 'symmetry' walls ws a-ia{P) = 13°" — (3'^~^ (a = 2, d) and 

^"The dominant linear forms can be identified in many pfiysically relevant cases with the simple roots of 
an hyperbolic Kac-Moody algebra [321 1331 [34] . 
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(ii) one curvature wall wid-idiP) = 2(3^ + jS"^ + ... + /J'^'^. Note that in d = 3, wus = 2/3^ 
corresponds to the (first of the) famous BKL walls of the form a'^ + 6'^ + c^, where 

2.3 Hamilton evolution equations 

Let us indicate, in a sketchy manner, the structure of the Hamilton evolution equations 
following from (jl5p . 

A 

where ca = ca{Q, P,dxP,d'^P,dQ,d'^Q), WAa denotes the (covariant) components of the 
linear forms WAiP) — WAaP'^- The system (jl7p is the one that we will analyze in detail in 
the sequel. 

Let us recall the basic classification of the set of dominant walls: either the fundamental 
chamber defined by the dominant inequalities wa{P) > is contained within the future 
(Sf=i f^"" > 0) light cone GabP°'P^ = 0, or it is not. The first case define what we call here 
chaotic systems, while the second defines non-chaotic systems. For instance, pure gravity 
in D = d + 1 is chaotic for d < 9 , and non-chaotic for d > 9. Note that this classification 
does not correspond to the often used asymptotically-velocity-terms-dominated (AVDT) 
systems versus non-AVTD ones. Indeed, there are AVTD systems that are chaotic. For 
instance, the Einstein-Maxwell system (in any dimension D = d + 1) is always chaotic and 
we shall see below that its asymptotic chaotic behavior can be described as a naive AVTD 
truncation of the full dynamics. Let us also mention that they are non-AVTD chaotic 
systems that are equivalent to AVTD chaotic system. For instance, gravity coupled to a 
{d — 2)-form is chaotic and non-AVTD because driven by its magnetic wall. However, by 
Hodge duality, it is equivalent to the Einstein-Maxwell system which is chaotic and AVDT. 

As a warm up towards understanding the structure of the evolution equations (jl7p we shall 
first consider the so-called non-chaotic systems. 



3 Iwasawa— variables treatment of non— chaotic systems 

In this section, we reformulate the results of \14\ I13| by using the Iwasawa variables, within 
an Hamiltonian approach. Let us recall that the treatment used in [141 [T3] consisted of 
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rewriting Einstein-matter systems into a Fuchsian form, i.e. 

drU — Au = e~^'^ f{x,T,u,dxu) , (18) 

where /i > and where the crucial conditions are (i) that the source term / should be 
bounded when r ^ oo (while the other variables take their values in a bounded set) and 
(ii) that the eigenvalues of the (space and time independent) matrix A be strictly larger 
than —n [35]; see appendix [B] for precise mathematical conditions. Then the main result of 
the Fuchs theorem is that there exists a unique solution u(r, x) of (jlSp which tends to zero 
as r — > 0. Moreover, the exponential decay of the source e~^'^ imposes a corresponding fast 
decay of solution which we shall write as 

= 0(e-^*"'") , 

where can be any number satisfying < /i*-^-* < fx [Note that can be as close as 
we want to 

Here we are going to show that the evolution equations in Iwasawa variables given by equa- 
tions ()17p can be rewritten in an alternative Fuchsian form which leads to a streamlined 
derivation of the results of [lll|13]. In order to do that, we need to do two things (i) define 
an asymptotic evolution system whose solutions {/3[o], ttjoj, Q[o], P[q]} parametrize the generic 
asymptotic exact solutions {(3,-k,Q, P}, (ii) rewrite the system of equations (fT7|) in terms 
of the differences u between {f3,TT,Q, P} and {f3[o],iT[o],Q[o], P[o]} such that the system of 
equations for u is Fuchsian, and (iii) define asymptotic constraints in such a way that the 
exact constraints are satisfied if the asymptotic constraints and the asymptotic equations 
of motion are fulfilled. This is done in the sequel and implies by the Fuchs theorem that 
there is a unique solution u that vanishes when r goes to infinity; the Fuchs theorem also 
tells us how u goes to zero as r ^ oo. This result gives a precise sense in which the ap- 
proximate solutions {/3[o], TTfo], Q[o], P[o]} parametrize the asymptotic behavior of the exact 
solutions {P, TT, Q, P}. If the asymptotic solutions {/3[o], vTfo], Q[o],P[o]} are general enough, it 
means that we have found the general asymptotic behavior of the gravitational field in the 
vicinity of a spacelike singularity. [More precisely, we want here that the solutions of the 
asymptotic system together with the associated asymptotic constraints contain the same 
number of arbitrary functions which is expected to enter the general solution of the exact 
constrained Einstein-matter system.] We will see that our new formulation is significantly 
simpler than that of [14:\ [T3] and is suggestive for approaching of the chaotic case. 

3.1 Definition of the asymptotic evolution equations 

The first step is to define a simplified system of equations that describe the asymptotic 
dynamics of the fields. There are, a priori, several choices for defining an asymptotic system 
when using Iwasawa variables. For instance, we can either neglect certain terms directly 
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= vr, 


drT^lO] a 


= 0, 


drQlo] 


= 0, 


drPu)] 


= 0. 



in the Hamiltonian or neglect some terms in the equations of motion. One of these choices 
gives a system essentiahy equivalent to the usually considered AVTD system in |14^ [13] . It 
would consist in keeping only the symmetry walls in the Hamiltonian (jlSp . Here, we will 
consider a technically simpler choice consisting in neglecting all the walls. Concretely, this 
means that we define the 'asymptotic Hamiltonian' as 

"^olAo]) Q[o]; TTjo], -P[o]] = ^G'^^7r[o]a7r[o]b , (19) 

the [(o)]'s refer to the zeroth order approximation of our general solution. The Hamilton 
equations corresponding to the Hamiltonian (jl9p are 



[0]i, ) 



(20) 

The solutions of these Hamilton equations are schematically (suppressing indices) , 

Ao] = PoT + Po , 
'''"[0] = Po , 

Qio] — Qo , 

Pm = Po, (21) 

where po, Po, Qo, Po do not depend on the time but depend on the spatial coordinates x*0 
Note that the metric corresponding to the 'asymptotic solution' (j21|) does not generically 
corresponds to a Kasner-type metric {i.e. a metric of the type ([T|)). Indeed, the Iwasawa 
'off-diagonal' variables AA's of a generic Kasner metric have limits as r — > oo but are r 
dependent for finite r, while the AA"j's corresponding to the solution (j2ip are constants [see 
paragraph 4.2 of [T^ for explicit expression of the Iwasawa variables of a Kasner metric]. 

3.2 Definition of tiie asymptotic constraints 

As the asymptotic Hamiltonian constraint, it is natural to take the asymptotic Hamiltonian 

W[o] = -G"^vr[o] aVr[o] b , (22) 

which has the useful property of being conserved along the asymptotic evolution equations 
(|2Up. Concerning the asymptotic momentum constraints, we need to know their structure 
in Iwasawa variables to be able to conclude. In view of this, we first express the momentum 



^'^As is usual when discussing Fuchsian theorems one makes the technical assumption that the spatial 
dependence of all the initial data {po{x), /3o{x), ...) is real analytic. 
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conjugate to the metric gij in terms of Iwasawa variables. 



Let T^i^a"''' := J\f°'iJ\f^j7r'^^ denote the Iwasawa-frame components of the momentum conju- 
gate to the Qij. Using this definition, the effect of the transformation of the configuration 
variables {gij} — > {/3",AA"j} on their conjugate momenta is obtained from writing 



ca 

iwa 



from which we can extract that (we recall that the metric in the Iwasawa-frame is g'^ 
e~'^^°'6ab so that 7ri„a"f, = e'^'^V; 



wa 

ab 

' ab\ 

wa / 

7i"iwa 6 = ~2'^b no sum over 0, 



^^cVTiwa'^a = ^'P\ only for i> a. 

In order to invert the above formula and get the TT-.^^'^a in terms of the tTq and V^a, let us 
rewrite the above equation for all i and a as follows, 

A/'""^*c(+)7riwa''a = 2'^'a[-] + ^(+) , (23) 

where X is a matrix defined by this equation, and where we have added to the various 
triangular matrices that appear an index referring to the fact that it is an upper/lower 
triangular matrix (+)/(—) or a strictly upper/lower triangular matrix [+]/[—]• We can now 
multiply equation ([23]) by = AA(_|_) and obtain, 

TT^Ja = i+)V\ [_] + M\ . (24) 

Let us decompose the matrix Hi^J'a into its strictly lower triangular part Tii^J'ay-]-, its di- 
agonal part TTi^J'b = —\'T^h and its strictly upper triangular part T^i^J'a[+\- The projection 
of both sides of equation (j24p on their strictly lower triangular parts yields an explicit 
expression for 7r''aiwa[-] (for h > a), namely 7r^aiwa[-] = \N'^i{^j^)V'^ a[-](^{b — a), where 

I 1 ifj;>0. 

Note now that TTi^J'ai+i, being obtained from the symmetric matrix vTi^j^"* = tti^J"^ by 
lowering an index by the metric g'iwaafc = e~^^"(5ab, can be related to iTi^J'a [-] in the following 
way, 

VTiwa a[+l = e ^TTiwa bl-] ■ (25) 



Finally, we have the following links 

b 



1 A 

vi"iwa b = — ^° overo. 



ifb>a VTi^a^a = TTiwa^'al-] = T^AA^i^^a , 



ifa>b = Oa[+] = ^e-2(/3"-/3^)AA^r,, (26) 
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Therefore, the T^i^^°'b[-\ are hnear in V^a and M'^a, while the vTi^^^fe [+] 's depend also on the 
/3's though the e"^^'^""'^''^ with a > b (see equation (I25p ). i.e. through symmetry walls 

Q-2wsbaW) ^ 







Let us now express the momentum constraints in Iwasawa variables 7ia. They read 

_l_ J_ c^i-.-ap •r-(p) 

where the Taai...ap^s are the F-'?' ^^'s expressed in the Iwasawa basis and the Tube's are the 
connection coefficients (with c denoting the differentiation index) in the Iwasawa basis, 

Note that vTi^^.^'^b is a tensorial density of weight 1. The Cj^a are the structure functions of 
the coframes Of^^, they are related to the structure functions C in the coframe u;* by the 
formula, 

- \a^.\c = dM\M-^\ - \M\C'jkM-^hM-^\ (27) 

Inserting these results in the expression for the momentum constraints gives 
1 1 

- = SfoTT a + C^^^'^cb'^ a + C'iwa acT^^ d " ^ (^"/^ )'^d 

+ i.^a,...a,_^(.) (28) 

Note that this is the general expression for the momentum constraints expressed in the 
Iwasawa variables. 



We then define the asymptotic momentum constraints by discarding the vTi^^ a [+1 contribu- 
tions in the exact constraint (j28p : 

~ '^jHayo] [db'^iwa a [— ] ~ '^daT^a + Ciwa cfe'^iwa a[—\ ~^ C'iwa ac TTiwa d[—] 

_^i_^ai...ap_^W^ ^jM^ no sum over a, sum over d (29) 

^^The comma in the expression gcd,h denotes the spatial derivative in the Iwasawa frame, i.e. gcd,b = 
fiiwa ij*ei(pcd) where Ci^^b = &i^a.b&i is the Iwasawa frame (dual to the coframe 9\^^: Siwa''(ei„aa) = <5a) and 
where Ci = et-'dj is the frame dual to the basic coframe = to^jdx^ used in equation Q. We will also 
sometimes denote gcd,b by dbQcd. 
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where the overall bracket [ means that one must do the following replacements vr — > 
7r[o], (3 — > P[o], Q ^ Qio], P ^ P[o]- Finally this definition corresponds (besides the replace- 
ment {(3,7r,Q,P} —>■ {/3[o], 7r[o], Q[o], -P[o]}) to setting to zero the symmetry walls e~'^^^"'~^''^ 
(with a > in the full momentum constraints. Along the solution (j2ip of the evolution 
equations (pOj) . the only time dependent term in Tiaio] is — 9a/?[J]7r[o] d/2 = —rdapt'^od/'^ so 
that we have the following relation, 

dr'Haio] = daTi-io] modulo equations ([20|) . (30) 

From this relation, we conclude that, when the Hamiltonian TC^q] = constraint is satisfied, 
the momentum constraints are conserved when the asymptotic evolution system (j20p is 
satisfied. Finally, it suffices to impose the constraints 7i[o] and TCaio] at any fixed moment 
to guarantee that they are satisfied for all time. 

Similarly, the asymptotic Gauss constraint for each p-form is defined to be the Gauss 
constraint with the asymptotic variables Q[o], -P[o] instead of Q, P: 

ai...ap_i — f) -TT "■i---<^P ^ n "1, -TT bo.2---ap ^ <^p-i, tt cii---bap 

'P(p)[0] ■— OapT^io] — 2'-'iwa[0] bapT^lO] — • • • — ^ '"'iwa [0] bapT^lO] 

+ Q^^^O]'"'apb7Tlor-"''''- (31) 

These constraints are preserved by the time evolution since the vr^]'^^"'"'' (which are some 
of the -P[o]'s) and the Ci„a[o]'^fec (which depend on the Af[o] G Q[o]'s via (i27|) ). are constants 
according to the asymptotic evolution equations (f20l) . 



3.3 Construction of a Fuchsian system for the 'differenced variables' 

Let us introduce the differences /3, vf, Q, P vi 

13 = Aoi+^, 
Q = Q [0] + Q 1 

P = Pioi+P, (32) 



^^where here /3[o], 7i"[o], Q[o] and P[o] are given in TTT^ . 
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and express the Hamilton equations in term of these variables. This gives 

drP--7t = 

+ a^(^e-2-A(/3[o])e-2«.A(^)) 
oP 

^ dQ ^ ddxQ 

_52(^e-2-^(^M)e--^(«) (33) 

where ca = ca(Q[o] + Q,Pio] + ^',5^(Ao] + P),dUPio] + ^),9(Q[o] + Q),d^Qio] + Q))- Let 
us sketch the proof that this system will be Fuchsian if all the 'walls' wa{P) entering the 
equation (jSSp are such that the following conditions hold, 

V^, Vx G U, wa{Po{x)) > e > , (34) 








1 

2 








\ 


































V 














/ 



where the po is the initial datum entering equation (j2ip which must also satisfy the con- 
straint (|22|l . i.e. C^^poaPob = 0, as well as the asymptotic momentum constraints ([25|) . 
Here, U denotes some open domain within the analytic d-dimensional manifold, on which 
one applies the Fuchs theorem. The fact that the system (jSSp is indeed of the form (jlSp for 
u = {P,Tt,Q,P) comes from two separated facts. First the matrix A being 



(35) 



is a nilpotent matrix and therefore is (thanks to the recent progress concerning Fuchsian 
systems |35) ) an allowed matrix A for a Fuchs system (see appendix [B]). Concerning the 
source term, let us show why the conditions (|34p guarantee that the 'source term /' — 
i.e. the right hand side of the system of equations (j33p — satisfies the right properties. 
Essentially the Fuchsian conditions boil down to requiring that the source term should be 
of order 0(e^^^) for some ^ > when {P, vf , Q, P, dxB, dx7t, dxQ, dxP, dxP, d^Q} take their 
values in a bounded set while r G [to, +oo] (see appendix |B]). The explicit time dependence 
of the source has three origins 

1. the Toda walls whose (r) time dependence is exponential e"^'"^(^M) = e-2"'A(po)T-2«;A(/3o) 

2. the various space derivatives appearing in the r.h.s. of (|33p can 'bring down', when 
operating on e~'^'^'^^P°^^^^'^ , one (dx) or two (9^) powers of r. 
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2. In addition, the dependence of the wah coefficients ca on SxAoi ^'^'^ ^xAo] means that, 
for some wallfl the coefficient CA can also involve one or two powers of r. 

Summarising, the r.h.s. of (j33p is a sum of terms of the form 

where P{t) is a polynomial in r. If we choose a /i strictly smaller than all the quan- 
tities 2ti;A(Po(a;)) > 2e > considered for any x £ U and any type of wall A, we can 
conclude that all the source terms in the equations (j33p are of the required order 0(e^''^) 
for r G [to, cxd] and x £ U. Therefore we can conclude that there exists a unique solution 
{(3{t, x),Tt{T, x),Q{t, x), P{t,x)} of equations ([55]) that vanishes when r — > cxd. Moreover, 
this unique solution satisfies (within the considered spatial domain U) the following estimate 
as r — > +oo: 



p = 


0{e 






TT = 


0{c 






Q = 


0{c 






p = 


0{c 




'1 



where < fi^ < /x being strictly smaller than the quantities 2wa{Po) > 2e). 
3.4 Constraints 

It remains to show that if the asymptotic equations of motions, and the asymptotic con- 
straints (P21[2UI[5T|) are satisfied, then the exact constraints ([TUIfTTlfT^ will also be satisfied. 

Let us first deal with the Gauss constraints (for notational simplicity, we shall consider the 
case of one p = 1-form, i.e. a Maxwell field). We recall, from equation (jl2p that, in this 
case, the exact Gauss constraint (in Iwasawa variables) reads, 

if := VaTT'' = daTT'' + Q^baTr" ^ . (36) 

This exact constraint is preserved by the exact equations of motions. We also recall that 
the asymptotic Gauss constraint ([3T|) . in our case, read, 

V[0] '■= VaTTfo]" = daTT[o]"' + Ciwa[ol''fea'^[Ol'^ ~ , (37) 

and is conserved modulo the asymptotic equations of motions (pO]) . We impose that the 
asymptotic Gauss constraint ([37|) hold. Then, the fact that the difference between the exact 
and asymptotic Gauss constraints are given by exponential 'walls' (entering the differences 
between the vr and ttjo] etc) implies that the exact constraint (|36p vanishes when r — > +00. 

^"^The walls depending on spatial derivatives of P are only 'subdominant' gravitational walls, see [12]. 
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Finally, from the obvious fact that an asymptotically vanishing quantity which is constant 
must be zero, we can conclude that the exact Gauss constraints are (weakly) fulfilled. In 
particular, the relationships V^T'^'' = are satisfied. 

Let us now consider the Hamiltonian and momentum constraints. To show that they are 
satisfied, we will argue that their evolution system is 'Fuchsian' in some generalized sense 
and therefore that there will be a unique solution of this system that vanishes when r — > oo. 
The Fuchsian system in question turns out to be homogeneous, so that the unique solution 
that vanishes when r — > +00 must be exactly zero. Then we will discuss why the exact 
constraints vanish when r — > +00 (when the asymptotic constraints are satisfied) and con- 
clude that, since they do, they must be the unique vanishing solution, i.e. zero. 



Let us first write down the evolution system satisfied by the exact constraints as a conse- 
quence of the Bianchi identities and the exact evolution equations. This reads 

d^n = e-2J^i'^'(v'^7^a-2^/?f„?^"), 

c 

drUa-Van = 0, (38) 

where the covariant spatial derivatives must take into account the weights of the various 
densities: 7ia has weight 1, and 7i has weight 2. The 'source term', i.e. the r.h.s. of the 
evolution equation for can be rewritten as, 

a c 

where Ha{P) = Ylbj^a ^ subdominant curvature wall (it corresponds to the special case 
of the curvature walls WabdP) with a = h). Therefore, the source term is an allowed one 
for the Fuchsian-like system (jlSp with u = {7i,7ia) and fi strictly smaller than all the 
quantities 2^a{Po{x)) > e > considered for any x £ U. However, this homogeneous sys- 
tem is not really a Fuchsian system because of the presence of the spatial derivatives term 
VaH in the second equation. The references [HI [13] tackled this problerr0 essentially by 
working with a suitably redefined Tl constraint, say Ti = e^'^TL. This redefinition produces 
for {'H,'Ha] a Fuchsian system if < < 2fia{po) (for all a and for all x £ U, Po{x) being 
as in equation (j21|) ). However, we think that there might be other ways of dealing with 
this problem. First, we expect that a generalization of the Fuchs theorem exists for linear 
homogeneous systems of the type psp . stating that the unique solution that vanishes when 
r +00 is everywhere zero. We then expect that one way to prove such a theorem is to 

^^See equations (|86|l and (|87|) in appendix [X] These equations were derived in a coordinate basis and 
without matter. It is obvious how to get the expression in a general basis. Moreover, we do not have to 
consider matter since V^T^^ = [as a consequence of the Gauss constraints (which are imposed to be 
satisfied) and the matter equations of motion]. Therefore, we can just replace in equations (|86|l and (|87p in 
appendix |X] the expressions (|80[) and (I81|l by their general expressions (|10|l and (lllf) . 

^® We translate here the argument used in [141 113] in our choice of variables and gauge. 
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work with an extended set of variables {Ti, Tia, Ga = dali.} and use the recent reference [35] . 
To summarise, the system (I38p is a Fuchsian-hke system that possess a unique solution 
(which is zero because the system is homogeneous) that vanishes when r +00. 

Moreover (i) as the definition of the asymptotic constraints 'H[o](/?[o] , 7r[o], (5[o], -P[o]) and 
'^a[oi(Aoi) ■^[0]) Q[oi) -P[oi) differs from the exact ones by neglecting some explicit exponen- 
tial walls in their mathematical expression, (ii) as /? — /?[o] etc tends to zero (essentially as 
0{e~^'^)) and finally (iii) as we have imposed the asymptotic constraints, we can conclude 
that the exact constraints 'H(/?, tt, Q, P) and TCaiP, ^) P) tend to zero as r — > +00. 

Finally, the constraints TC and Tia being uniquely defined as being an asymptotically van- 
ishing solution of an homogeneous Fuchsian system, vanish for all times r. 

Summary: An asymptotic solution {/?[o] (r, x), ttjo] (r, x), (^[o] (t, a;), P[o] (r, x)} ([HP obeying 
the asymptotic constraints [2^1 . the asymptotic evolution system (2D(l and conditions (3^ , 
parametrizes a solution {/3(r, x), 7r(r, x), Q(t, x), P(t, x)} of the full constrained 

Einstein-matter equations (this is pictured in Figure\^. Moreover, the asymptotic closeness 
of the two solutions satisfies inequalities of the type, 



P{t,x] 


1 = /? 




= 0(e" 






Q{t, x) 


= Q- 


~ Qio] 


= 0(e" 




-)r 


7t{T,x] 


1 = TT 




= 0(e" 


-/.{- 


-)t 


P{t,x) 


= P 


- -P[oi 


= 0(e- 


V- 





(39) 

where fi is any number strictly smallerMx G U andWA than the quantities 2wAiPoix)) > 2e, 
when T +00 . 

On the generality of the construction 

An important fact to notice about the above construction is that the solution of the asymp- 
totic equations ([2T]) subject to the constraints ([22|) . ([29]) and asymptotic Gauss constraints 
possesses as many arbitrary functions as one expects to be present in the 'general solution' 
of the constrained Einstein-matter equations. The inequalities ([Mj) impose restrictions on 
these arbitrary functions but do not change their number. Let us repeat that this con- 
struction applies only in cases where the fundamental chamber in /3-space defined by the 
inequalities WAiP) > extends beyond the light cone Gab/3'' < 0, /? > as illustrated 
in Figure [2j 
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Figure 2: Non-chaotic behavior This picture is a schematic drawing of the asymptotic dynamics 
of the 'diagonal variables' at a given spatial point x, this dynamics is represented in the /3~space. 
The 'dashed' arrow represents the asymptotic solution /3[o] (which is valid after the last collision 
on a wall and corresponds to a free motion of the particle /3). The exact solution is sketched as 
a continuous curve. The idea is that the approximate solution /3[o] becomes better and better as 
T — ^ +00, this is formalised by the Fuchs theorem that tells us precisely how f3 — /3[o] when 
T — > +00, see the text. Note that here we consider a non-chaotic system and that the 'fundamental 
chamber' determined by the walls in not contained within the light cone. 




Comparison with the asymptotically velocity dominated system 

The appendix [Cl compares the velocity dominated system of |141 I13] and our approach. The 
essential differences between the two approaches are the following: (i) they do not use the 
same asymptotic system, (ii) the use of Iwasawa variables allows for a more transparent 
treatment of the 'source terms' in the Fuchsian equations (indeed, in Iwasawa variables, it 
suffices to read off the exponential terms in equations (j33p ). (iii) the use of Iwasawa vari- 
ables avoids the technical problems linked to measuring the 'difference' between the exact 
and asymptotic metrics by writing {g^^g)"'b = ^""b + t""''7"b with some carefully chosen a's 
and an asymmetric matrix 7. 

We should, however, remark that the two methods differ in the extension of the open 
regions U where the Fuchsian method can be used to construct the metric. In the method 
of |14l [T3] one can cover the full analytic manifold by using many small neighborhoods in 
which the frame approximately (but analytically) diagonalises the second fundamental form 
K. In the Iwasawa approach one can work in large open domains, but there is a problem 
connected with the presence of co-dimension 2 submanifolds where two eigenvalues of K 
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coincide. This problem is briefly discussed in appendix [Cl More work is needed to extend 
the Iwasawa-variable approach so as to be able to cover the full analytic manifold. 

4 Iwasawa— variables treatment of chaotic systems 

We now turn to our main purpose, which is to give a precise formulation of the asymptotic 
BKL behavior in the chaotic case. In this perspective, we follow the same strategy as for 
the non-chaotic case: 

(i) We first define an asymptotic evolution system, made of ordinary differential equations 
(ODEs) that will describe the generic asymptotic unconstrained dynamics of the Iwa- 
sawa variables near a spacelike singularity. Of course, the solutions {/9[o], ttjoj, Q[o], -P[o]} 
of the asymptotic system are much more involved than a Kasner-like behavior and 
cannot be given in a closed form. 

(ii) We then define asymptotic constraints whose vanishing is preserved by the above 
defined asymptotic evolution system. 

(iii) Next, we construct a ^generalized Fuchsian system'' that describes the behavior of the 
differences (3 = P — /3[o], vf = vr — ttjoj, Q = Q — Q[q] and P = P — P[o] between an 
exact solution {(3, vr, Q, P} of the considered Einstein-matter system, and a solution 
{/9[o], TTfo], Q[o], -P[o]} of the asymptotic evolution system. We then argue that, given a 
solution {/9[o], TTjo], Q[o], -Pp]} of the asymptotic evolution system, there exists a unique 
solution {(3, vf , Q, P} of the differenced system which goes to zero as r — > +00. 

(iv) We formally show that, if the asymptotic constraints are satisfied, the full constraints 
satisfy a generalized Fuchs system. We then argue that the full constraints will be 
satisfied as a consequence of the vanishing of the asymptotic ones. 

Finally, our methodology suggests that one can indeed parametrize a solution of the full 
constrained Einstein-matter system by a solution of the system of ODEs defined in (i). 

Note: we use the same notation for the solution of the asymptotic system of equations, 
asymptotic solutions, asymptotic Hamiltonian etc in the chaotic and non-chaotic cases. 

4.1 Definition of the asymptotic evolution equations 

The billiard picture provides a guide for choosing a suitable asymptotic evolution system 
since it gives us an intuitive description of the asymptotic dynamics. In the billiard approxi- 
mation, the dynamics of the 'diagonal variables' /3's is described as a free Lorentzian motion 
interrupted by refiections upon infinite-potential walls and the 'off-diagonal variables' Q's 
are frozen. Here, we shall go beyond this simplified 'sharp wall' billiard picture and work 
with exponential ('Toda') potential walls. The asymptotic system should completely deter- 
mine the asymptotic dynamics, given some suitable initial data. It is therefore crucial to 
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use a system of ordinary differential equations (rather than partial differential equations) to 
characterize the asymptotic dynamics, so that we are able to use theorems about existence 
and uniqueness of solutions. 

Motivated by these reasons, we define (for any chaotic Einstein-matter system) an asymp- 
totic evolution system in the following way: for the 'diagonal variables', we keep in the 
Einstein-matter equations (|17p only the dominant exponential walls while for the 'off- 
diagonal variables', we neglect all the walls in the equations of motion. These prescriptions 
define an 'asymptotic evolution system' which reads, in sketchy form: 

A 

drQlO] = 

drPm = 0. (40) 

Here A labels the dominant walls only. As exhibited in equation (jl6p the coefficients c^i of 
the dominant Toda walls depend only on the Q's, P's and the first spatial derivative of the 
Q's. It is crucial that they do not involve the spatial derivatives of the /?'s. The solutions 
of the last two equations in the system ([10}) are simply (5[o] = Q(^o){x), P[o] = P(o)(2;). 
Considering Q(q){x) and P(o)ix) as given data, and replacing them in the other equations of 
the system PU]) . we see that the diagonal variables {/3[o], ttjo]}, at each given spatial point x, 
satisfy a system of ODEs. It is then easily checked that the latter system of ODEs follows 
from the Hamiltonian Tio'. 

^oQ(o),P(o)][Aoi;^[oi] = ^G"^'vr[o,aVr[o,6 + 5]c^(^'(o),Q(o),5aQ(o))e-2--^(^m) . (41) 

A 

Note that, from [12], the qualitative behavior (as r — > +oo) of the solution for such systems 
of equations is as follows: the vr^i's go to zerc^ZI and the /3[o]'s behave approximately as in 
the sharp billiard picture (free motions 'pr + consf interrupted by collisions against the 
'walls' e-^"'^^^)). 

4.2 Definition of the asymptotic constraints 

It is natural to define the asymptotic Hamiltonian constraint to be 

W[oi := ^G'^'vrioiaVTioift + ^c^(^'[o],Q[oi,5aQ[o])e-2-^(^m) . (42) 

A 

Like in the non-chaotic case the asymptotic Hamiltonian constraint (j42p coincides with the 
asymptotic evolution Hamiltonian ()4ip . and is therefore preserved by the asymptotic time 
evolution. Let us now define the asymptotic momentum constraints as the formal limit of the 

Actually, this property is guaranteed only if one imposes the constraint 'Ho Q,g,p,gs ~ 0. 
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'full' momentum constraints when r +00. We start from the general expression (|28p for 
the momentum constraints. In view of what was recalled from [12], vr — > as r ^ +00 (after 
imposing (j42p ). so that we can discard the terms linear in vr in the previously considered 
non-chaotic asymptotic momentum constraints ()29p (which neglected exponential walls 
that we still formally neglect). We therefore define the asymptotic momentum constraints 
as follows, 

~ 2'^0,lO] ■~ [db'^iwa a [ — ] + C'iwa cfe'^iwa a [— ] + Ciwa aC^Tiwa d[ — ] 

_^i_^ai...ap_^^p)^ ^j[0]^ no sum over a, sum over d (43) 

where the J-'aai...ap^s are the -?^{'5|i^)...(ip)'s expressed in the Iwasawa basis, where vTfoj^a [-] is de- 
fined as the r.h.s. of the second equation (I26p and where the overall bracket means that 
one must do the replacements Q Qio], P Pio]- Note that the momentum constraint (j43p 
contains only the time-independent quantities A/jo] , S^A/jo] , Q[o], P[o] and therefore is trivially 
preserved by the time evolution. 



Finally, the asymptotic Gauss constraint for each p-form is defined to be the Gauss con- 
straint with the asymptotic variables Q[o], P[o] instead of Q, P, i.e. 

ai...ap-i ._ Q ai...ap ^ ai ba2...ap ^ f~i a.p-1, „ ai...bap 

f(p)[0] ■— C'ap^m — 2^iwa[0] bapT^lO] — ■•• — -l^iwa[o] bupT^lO] 

+Q^m"''apb7rior (44) 

These constraints are preserved by the asymptotic time evolution since the 7r[o]"^'""'' (which 
is one of the -P[oi's) and Ci^^io]"'bc (which depends on the M[o] £ Qpi's via (f27|) ). are constants 
according to the asymptotic evolution equations (POj) . 



4.3 Construction of a 'generalized Fuchsian' system for the 'differenced 
variables' 

We now rewrite equations (jl7p in terms of the differences {/?, vf, Q, P}, 

T^a = '''"[o] a-\~ T^a 1 

Q = Qio] + Q 1 

P = Pm+P. (45) 
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This gives (when suppressing indices on (3 and vr), 
9./3-^7f = 

A 
A' 

+ Qj^^e-'^A'iPm)f,-'i^wCP)) 

A' 



dP 
A 



-2«,_4, (/3[o] ) -2w^, {^) 

A- 



dQ 'dd^Q 
dcA' 



_ ^a2(^e-2-^'(Aoi)e-2-^'(/5)) , (46) 

where we recall that A labels the dominant exponential walls (with coefficient c^^ which 
depend only on {P, Q, daQ}) while A' labels the subdominant exponential walls (with coef- 
ficient c^/ which depend on {P,Q,daQ,d'^Q,daP,d^P}). In addition, in all the coefficients 
and c^/ on the r.h.s. of the system (j46p one must do the replacements (j45p . so that 

for us CA = ca{Q[o] + Q,P[o] + P,d{Qio] + Q)) and ca' = c^'(Q[o] + Q,Pio] + P,dxiPio] + 
92(/?[o, + 9(Q[„j + Q), a2(Q[o] + Q)). 



This system of equations is not a Fuchsian system as defined above. However, it is similar 
to such a system. Indeed, it contains a space and time independent matrix A which is 
again given by (|35|) . However, the crucial difference between (|46|) and a Fuchsian system 
concerns the source term on the r.h.s.. Instead of containing (modulo a bounded term) a 
space-independent factor which is exponentially decreasing with r, e~^'^ , it contains expo- 
nential wall terms e"^^-^^'^!"!-' and e"^"'-^' ^'^["1 where /3[o](t, x) is a solution of the asymptotic 
evolution system (P0|) . 

Let us qualitatively analyze the behavior of this source term as r ^ +00. For that, let 
us start by recalling the sketchy time dependence of the solution /?[o] that would be given 
by the billiard picture: namely a succession of Kasner epochs. During each Kasner epoch 
the source term is exponential decreasing (in the time coordinate r), indeed, during the 
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'Kasner free motion' /3[o] ~ PoT + Po so that each exponential wah factor term e~^"'^*-'^[°l^ 
{A = A, A') behaves has e~^^^^'^ with fi{x) ~ 2wa{Po{x)). However, this exponential de- 
crease is interrupted around the instants of collision on the dominant walls, during which 
w_a{Pio]) in fact vanishes so that e"^"'"*^^!"!-' would seem to become unity. By contrast the 
subdominant terms e"^'"-^'^'^!''!-* are always exponentially decreasing because the /3-particle 
generically never hits them, being deflected by a collision on dominant walls before reach- 
ing them. [We are here neglecting the measure zero set of trajectories which exactly hit a 
'corner' of the billiard, where a subdominant wall intersects dominant ones.] 

The sharp billiard picture just recalled is only an approximation to the asymptotic dynam- 
ics (|40p . When taking into account the existence of exponential walls in the equations (j40p 
one can describe more precisely the behavior of (/?[o], vr[o]) and thereby of the crucial 'source 
terms' oc e"^'"'*^^!"'-* appearing on the r.h.s. of (j46p . Indeed, following the method used in 
[12| one can conveniently analyze the dynamics of (/?[o] , vr[o] ) following from the asymptotic 
Hamiltonian Tiio] (j^Tj) . and submitted to the zero-energy constraint (P2]) . When decompos- 
ing (3°' as (3°' = with = —GahP°'(i^ and Gahl"'!^ = — 1, which correspondingly implies 
tTq = TTa — vrp7(j where lil (submitted to the constraint ^""nZ = ) is conjugate to the 
'position' 7*^ on the unit hyperboloid {Gahl""!^ = and where VTp is the conjugate to the 
variable p, one finds that the Hamiltonian reads 

W[oi = ^(-vr^ + ^) + V[oi, (47) 

where 

V[o,(/9,7):=E^-4«""'""^'^- (48) 
The zero-energy constraint can then be written as 

i(-7r2+7r2)+p%=0, (49) 

where tta := p'Kp is now conjugate to A := In p. From this constraint one infers (see |12j ) 
that, as r ^ +oo and therefore p — > +oo, vr^ tends to a finite limit say and therefore \tt^\ 
oscillates between px (far from the walls) and (during a 'collision'). From this result one 
also infers that the maximum value of Wfo], reached during a 'collision' (i.e. when |7r^| = 0), 
is such that p^Vyo] = \t^\ \p\ as r — > +oo. Finally, one concludes that as r ^ +oo, and 
therefore p — > +oo (roughly proportionally to r) even the maximum values of the domi- 
nant exponential potential (reached during the collision) decay like V[o] oc p^^. [One also 
conclude from tTq = p^^ (vrj — 7rA7a) that the components of the /3-conjugate momenta vTa 
decay proportionally to p^^-] 

Summarising, we conclude that each dominant potential term e~^"'-^(^[oi) entering the r.h.s. 
of (j46p has the qualitative behavior depicted in Figure [3l namely an overall exponential 
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decay, interrupted by 'peaks' (of decreasing magnitude oc corresponding to collisions. 
In addition, (see appendix A of [l2]) the r-time spacing between successive peaks increases 
(roughly like In/) ~ Inr) as r — > +00. As for the behavior of the subdominant exponential 
potential terms e~^"''4'^^M-' entering the r.h.s. of (06]) it is expected to be somewhat similar 
to the one depicted in Figure O except for the facts that the overall exponential decay 
should be faster, and that the peaks should be much rarer (corresponding to a collision 
happening nearly in a 'corner'). [We expect that the faster decay and the rarer occurrence 
of peaks also compensates the fact that the presence of {dxjS)"^ and 9^/3 in the coefficient 
c^' generates a growing behavior oc of the c^'.] 

fir) * 




Figure 3: Schematic drawing of the source term of the system of equations (|46)) . 

The system (flSj) for {^, 7f,Q,P} (in which the r.h.s. depends on a solution {/?[o], VTfo], 
Q[o])-P[o]} of the asymptotic system (|40p ). can be viewed as a generalized Fuchs system. 
Note that the structure of this 'generalized Fuchs system' is of the form 

dru -Au = Y, e-'"-^(''m Va(Aoi, ^xAoi, ^, d^u, dlu) , (50) 

A 

where m is a vector-valued unknown function u{t,x) = (Mi(r, x), Mjt(r, x)), the linear 
forms wa{(3) are the same ones that enter in the system (|46p and where the source terms 
/a can be read off the system ()46p . In view of the arguments given in [12] (in particular, we 
recall that in the appendix A of this reference, it has been argued that the 'peaks' in the 
source terms pictured in Figure [3] are such that their integrated effect allows u to have a 
limit as T ^ +00) and partially recalled above, we expect that (under the conditions spec- 
ified below) there exists a unique solution {P, n, Q, P} of (06]) tending to zero as r — > +00 
(and more generally a unique solution u of the system (|50p . given suitable conditions on /, 
which tends to zero as r — > +00). 

The conditions necessary for this result to hold for (j46p are expected to be the following: 

• the asymptotic initial data Qo{x), Po{x) must be such that the coefficients c_a{Qo{x), 
Po{x),dxQo{x)) of the dominant potential walls remain strictly positive over the con- 
sidered domain U, 
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• the asymptotic initial data Po{x) = /3(ri,x), tTo^x) = tt(ti,x) at some finite time ri 
must satisfy the zero-energy constraint 7i[o]i(3o,'^o,Qo, Po) = 0. 

4.4 Constraints 

To complete the story, we need to check that the exact constraints are satisfied (along the 
exact equations of motion) once the asymptotic ones are fulfilled (along the asymptotic 
equations of motion). The reasoning is the same as in the non-chaotic case: 

1. First we treat the Gauss constraints. We impose that the asymptotic Gauss con- 
straints (j44p hold. It is again obvious that these asymptotic constraints are preserved 
by the asymptotic equations of motion (j40p . On the other hand, the exact Gauss 
constraints (|12p are preserved by the exact evolution equations. Moreover, they differ 
from the asymptotic ones by exponential 'walls', and consequently they vanish when 
r +00. As in the non-chaotic case, we can conclude that they are equal to zero 
(because they are constant and they tend to zero). 

2. We then turn to the Hamiltonian and momentum constraints. Let us require that 
the initial data Qo, Po of the asymptotic evolution system satisfy the asymptotic 
momentum constraints Ti.[o]a = 0, (j43p in addition to the asymptotic Hamiltonian 

(■[0] 



constraint "Hrgi = 0. The constraints obey the evolution system ([38|) F^ i.e. 



drU = j;e-2^«(^)(9,W, + 2/3'^,W,-2(^/3^),,H,), 

a c 
drHa-VaH = 0. (51) 

As in the non-chaotic case, the above system is not Fuchsian due to (i) the term Vq'H 
[however, as we have argued previously, this term should not be a problem], (ii) the 
source term (i.e. the r.h.s. of the first equation in the system (jSip l which is not an 
allowed Fuchsian source term (this fact contrasts with the non-chaotic case) . However 
we have already dealt with this kind of source term in the system (|46p . In the present 
case, the source term is even 'better' than the one of (j46p because it contains only 
'subdominant walls' which decay faster than the dominant ones and which exhibit 
rearer 'peaks' (because the 'peaks' occur when the 'ball' hits a corner) [however these 
'peaks' contain a factor oc r^] (see discussion in section . Accordingly, the system 
(jSip is a 'generalized Fuchsian' system in the sense given in section 14.31 and we con- 
sider it likely that it possesses a unique solution that vanishes when r — > -|-oo. On 
top of that, the system (j5ip is homogeneous and this implies that the unique solution 
in question is zero. 

We can then conclude that the exact Hamiltonian and momentum constraints are 
satisfied since (i) they differ from the asymptotic ones, which are imposed to hold. 



^See appendix |X] and the discussion in the section about the constraints in the non-chaotic case 
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by exponential 'walls' (this implies that they decay as r — > +00), (ii) they obey a 
'generalized Fuchs system', as just argued, and therefore there exists a unique solution 
that vanishes when r — > +cxd. 

Summary: Let us summarise our conjectural results concerning the asymptotic dynamics 
of the fields in the vicinity of a spacelike singularity for a chaotic Einstein-matter systems 
in the following statement: 

Let 

• {Q[o]{^)i be functions of the spatial coordinates such that the coefficients 
ca{Q{o]{x), P[o]{x), daQ[o]{x)) nowhere vanish, and that the 'fundamental chamber' de- 
fined by the inequalities w_a^{(3) > is contained within the future lightcone GabP°'P^ = 

• (/3[o] , TTfo] ) be a solution a the asymptotic system of equations ([^01) with initial condi- 
tions (/3[o](x),7r[o](a;)) given at some finite r = ri which satisfy the asymptotic Hamil- 
tonian constraint (j42|) . and 

• impose that the asymptotic momentum constraints (|43|) are satisfied at r = ri as well 
as the asymptotic Gauss constraints (HH). 

Then there exists a unique solution (/J, vr, Q, P) of the Iwasawa-variable form of the full 
constrained Einstein-matter equations such that the differences (3 = [5 — /3[o] tt = tt — vTfo] Q = 
Q ~ Q\o]iP = P — P[o]} tend to zero as r — > +00. 

5 Pure gravity in dimensions 4 < D < 10 

To give a more concrete example of the general formulation of the BKL behavior of Einstein- 
matter systems discussed here, let us consider the specific example of pure gravity. We 
consider spacetime dimensions D such that 4 < D < 10, so that the corresponding behavior 
is generically chaotic. For this case, our precise formulation of the BKL conjecture is the 
following (we denote hy d := D — \ the space dimension 3 < d < 9): 

(i) Initial data Let us give ourselves the following initial data: 
d{d — l)/2 spatial functions J\f(o)"'i{x) for a < i, 

d{d — l)/2 spatial functions P(o)*a(a;) for a < i, 
d spatial functions /3(o)"(x) and 
d spatial functions TT^o)a{x). 

(ii) Asymptotic Hamiltonian Given these data we define the following asymptotic 
Hamiltonian 

W-y-P(Aoi,vr[o,) = ^G"^'7r[o]a7r[o,6 + Vr°^P + V^^^'^P, (52) 
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Figure 4: Chaotic behavior This picture is a schematic drawing of the asymptotic dynamics of 
the 'diagonal variables' at a given spatial point x, this dynamics is represented in the /3-space. The 
dashed curve represents the zeroth order solution /3[o] . The exact solution is sketched as a continuous 
curve. The idea is that the approximate solution /3[o] becomes better and better as r ^ +00, this 
is formalised via a 'generalized Fuchs theorem', see the text. Note that here we consider a chaotic 
system and that the 'fundamental chamber' determined by the walls in contained within the light 
cone. 




a=l 



(where i = 1, . . . , d is summed over) and 

VG^'^P = ^e-2— (/^[o,)(qi,^_,,)^ (54) 

In the last equation, OabciP) (for b ^ c) denotes the linear form OabciP) = + 
13^ (evaluated for a = 1, 6 = d — 1 and c = d), and C(o)"fec (with h ^ c and 

C{o)"f)c = — C(o)"cfe) denote the structure functions {dO^Q)"- = — ^ C(o)"6c A 0(o)^ ^(o)'^) of 
the 'asymptotic Iwasawa frame' 9^o^^{x) = M{o)i"'{x)u)^ ■ Note that all the coefficients 
entering the exponential potential terms ([53]) and depend only on the spatial point 
(through 'P(o)(x), M(o){x) and dxM^o){x) which enters C(o)), so that the asymptotic 
evolution system for /? and vr constitutes, at each point of space, a well-defined system 
of ODE'S. 

(iii) Asymptotic evolution equations The equations of motion deduced from the Hamil- 
tonian (j52p are called the asymptotic evolution equations, they are the 'chaotic analog' 
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of the AVTD evolution system considered in the non-chaotic, monotonic power-law 
case. They are of the form: 

d 

O P[o] 

+ V™(Aoi;7'(o),AA(o),9.AA(o))] . (55) 

(iv) Asymptotic constraints We impose that the initial data satisfy the following asymp- 
totic constraints, 

w^^y'^p ( Ao] , vTio, , AA(o) , AA(o) , P(o) ) = , 

Wr-P(AA(o),a,AA(o),P(o)) = 0, (56) 

where 'H^^'^'^ is the (conserved) quantity defined in equation ()52p . and where the 
definition of TiT^™"^ is equation ([33]) above. 

Finally, this leads to the following precise formulation of the: 

BKL conjecture in Iwasawa variables. Let, for x £ U, the spatial functions V^o) -^(o) 
and C(o) be such that the d x~dependent coefficients V{o)^aJ^{o)"'~^^i and C(o)^d-id (whose 
squares define the coefficients of the d exponential potential terms ([53]) and ([M]) ) do not 
vanish in U. Let (Aoj (''"i 2;), 7r[o] (r, x), be the unique solution of the asymptotic evolution 
system ([55]) with initial conditions Aoi(''"ii2^) = Ao)(^) and 7r[o](Ti,x) = TT(o)ix) at some 
finite time r = ri and satisfying the asymptotic constraints ([56p . Then there exists a unique 
solution {P{t, x),-k{t, x),M{t,x), V{t,x)) of the vacuum Einstein equations (including the 
constraints) such that the differences (3{t,x) = f3{T,x) — (3[o]{t,x), 7f(r, x) = 7r(T, x) — 
7r[o](r, x), M{t,x) =M{t,x) — AA(o)(x), V{t,x) = V{t,x) — 'P(o)(x) tend to zero as x £ U is 
fixed and r +00. 

6 Asymptotic geometrical structure on cosmological singu- 
larities 

In the previous sections, we studied in detail the asymptotic dynamics of the gravitational 
field in the vicinity of a spacelike singularity in Iwasawa variables. In particular, we have 
seen that some of the variables have limits when r +00, i.e. the Q's and P's while the 
/?'s have no limits. Moreover, we have argued that in the 'chaotic' case the Tr^'s tend chaot- 
ically to zero. Of course these Iwasawa variables are dependent on the choice of coframe 
(J* := uj^j{x)dx^ used in the equation ([3]) [where a;* could be simply a coordinate coframe]. 
This raises the question to know whether the Iwasawa variables, despite their 'gauge depen- 
dence' capture some well defined geometrical structure at the BKL limit and what is this 
structure. In the non-chaotic case, this question has a clear answer. Indeed, the frames that 
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diagonalise the second fundamental form with respect to the metric have a well defined limit 
at the singularity. They therefore provide a fields of 'directional frames' at the singularity, 
i.e. a field of frames considered modulo rescalings of each frame vector. In this section, we 
will investigate the chaotic case. 

Our starting point is the existence of many variables having finite limits at the singularity, 
namely the AA^j's and V\'s, (say Af^iir, x) M{o)°'i, {x) V^ir, x) V(o)\{x) as r ^ +00). 
The problem is, however, that the quantities M(o)"'i and 'P{o)*a do not have, a priori, a clear 
geometrical meaning because they depend on the coframe used on M^. One way of 
addressing this issue is to act on the coframe (j* with an arbitrary local transformation 
A G GL{d,M) to investigate what information is left invariant in A/'(o)°j(x) and V(o)^aix) 
when 'rotating' the coframe by an arbitrary A. We will try to assign canonical values to 
the Iwasawa variables M(o)"'i{x) and V(^o)'a{x) at the singularitj^ by means of a suitable A 
(e.g. in the non-chaotic case we can find A's such that: J\f(o)°'i{x) and V{o)^a{x) 0). 

If we are able to assign canonical values to all asymptotic values of the M(o)'s and 7^(o)'s, 
this would mean that we would have again privileged directions at the singularity like in 
the non-chaotic case. We will see that the situation is actually more subtle than this. 

A general matrix A € GL{d,M) can be decomposed into three parts (i) a diagonal part, (ii) 
an upper triangular matrix (with ones on the diagonal) matrix and (iii) a lower diagonal 
matrix (with ones on the diagonal). 

(i) The action of the diagonal part of A consists in shifting the values of /5's. However, as 
the /3's have no limit as r ^ +00 it is not clear how to extract some geometrical meaning 
from such shifts of the /3's. 

(ii) Concerning the action of the upper triangular part, let us show that it can be used to 
fix the asymptotic values of the ^^{oj^j's to be 5^i. The crucial point is that if we act on the 
coframe by an upper triangular matrix, i.e. lo'^ = A^jlo^ , since 

a 

must be equal to ds''^ = e"^^"'AA"'iAA'" iuj'^uj'^ and since the transformation is defined 
by demanding that N and N' be both upper triangular, we easily see that 

M"'^ = M"^jA^i (for upper triangular A only) . 

This result is actually valid for any r, and yields, in the limit r — > +00, J\f^o) = AA'o)A. 
Now it suffices to perform the transformation with A*j chosen to be the upper triangular 

'^^Note that, in the non-chaotic case, we have access to more information, namely the non-zero limits of 
the tt's. 
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matrix A* 



A/(o)'j to fix the canonical values of the AA'^j's to be the unit matrix. Note 



that after this fixing of the frame (such that AA'^j = 6) the frame a;'* becomes identical at 



the singularity with the limiting Iwasawa frame 9'^ = AA'". 



UJ 



UJ 



(iii) Let us now consider the effect of the remaining freedom in a general A G GL{d, 
i.e. a lower triangular matrix, namely, 



A 



/ 1 

A^i 1 












1 

A^d-1 








The action of any A (upper or lower triangular) on the frame components of the metric 
ds^ = gijio^Lo^ = g^jUj'^Lo'^ is always given by the following linear action, 

<7^(r) = A-i^5fcKr)A-^- . (57) 

However, the induced action of such a A on the Iwasawa variables /?(t) and AA(t) parametriz- 
ing gijir) is somewhat complicated when A is lower diagonal because it is non linear (con- 
trarily to the simpler case just discussed, of an upper triangular matrix whose action was 
nicely compatible with the upper triangular nature of M and was thereby linear). Let us 
then consider the case where the remaining lower triangular matrix A is close to the identity, 
say 

A = 1 + A, 

with infinitesimal (strictly lower triangular) A. In addition, as we always assume that we 
have already used an upper triangular matrix A+ to fix A/(o) (after dropping the primes) to 
the identity, we can write that AA(t) is of the form. 



Mir) 



l + n(T), 



with n(r) ^ as r — > +oo. As both matrices A and M{t) can be treated as infinitesimal 
elements, we will neglect terms of order 0{n'^) and O(A^). Note also that, as we are near a 
situation where Af^i ~ (5?, the distinction between the 'a' type indices and 'i' type indices 
disappear. When replacing the Iwasawa decomposition gij = e~'^^°'J\f"'iJ\f"'j in equation 
(j57p . we find at the linear approximation, 

= n'j + [n,XYj + Oie-'^^'-^'y), 



n 



-2/3" 



e-2/3'(l + 2[n,A]^ 



no sum over i . 



(58) 



As the P's are the canonically conjugate to the N = 1 + n, the law of transformation of the 
P's is obtained by using the 'conservation' of the canonical form VdM+irdP = V dM' +tt' d(3' . 
We then easily find 



V'j + [Wfj - X'jW^ - vr'*) + 0{e 
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In addition we see from equation (j58p above that the hmit n(r) — > is invariant under 
such a lower triangular transformation A, so that the canonical value M(o) = ^^^^ fixed 
by such a A. Since at the BKL limit, the exponential 'symmetry walls' e~^*^^^^^') (j > i) 
vanish and that, in the chaotic case, the momentum conjugate to the /3's is going to zero 
(after each collision, they are 'redshifted' [E]), we obtain the following action of the lower 
triangle A = 1 + A on the limiting values of the Vs, i.e. 7^(o)'s: 



' (0) ] 



^(o)^- + [A,P, 



(0)J J ■ 



(59) 



In view of the strictly lower triangular nature of both A and V it is easily seen that the 
transformation law (|59|) implies that the elements of V(^o) on the first lower diagonal are left 
invariant: 



' (0) » 



no) I 



Consistently with our general requirement of having non vanishing dominant symmetry wall 
coefficients, we assume that all the constants V(o)^~^^i are non vanishing^ Then it is easily 
seen that by choosing a suitable A one can change at will the values of the 'P(o)'s on the 
lower diagonals, 'P(o) {n > 2). This proves that there exists a A such that we can fix 
P(o) to the following canonical form. 



Vo 



/ 














\ 




v\ 


















r\ 


























V 








... V'd-i 





/ 



(60) 



Let us now study what are the coframe changes that leaves the canonical of M^o) (AA(o) = S) 
and P(o) ('P(o) = 'Po) invariant. We already know that the canonical form J\f^o) = ^ fixes the 



upper part of A to be unity. As for the lower diagonal part A 
be fixed to Vo implies, from the equation ([59]). the condition 



1 + A, the request that 7^, 



(0) 



0. 



(61) 



To analyze the consequences of this constraint, we decompose A into a sum of matrices with 
non vanishing elements only on one of its 'lower' diagonals, i.e. A = Ai + A2 + ••• + A^^^i 
where the only non zero elements of A„ are (A„)*^"j. It is then easily seen that Ai must be 
proportional to Vo- Then one similarly finds that A2 oc Vo etc... Finally the most general 
A = 1 + A fixing 'P(o) to its canonical value Vo is found to be of the form. 



A = aiVo + 02^0 + asT^o + - + "d-i^, 



d-l 



(62) 



for some constants a^. 



Such a nilpotent element P(o) is called 'regular' in the mathematical literature. 
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Let us discuss the geometrical meaning of our findings. We consider again a general A G 
GL{d, M), containing both upper and lower triangular parts (for the reasons explained above 
we do not consider the diagonal part). If we had been able to define a canonical form whose 
'stabiliser' in GL{d,M) had been only the unit matrix, this would have meant the existence 
of a preferred directional frame (frame modulo rescalings) at each spatial point x 'on the 
singularity'. On the other hand, if the stabilizer had been the full group GL{d,M), this 
would have meant that no preferred frame at all remained on the singularity. Actually we 
have a intermediate situation, our stabiliser 5" is a proper subgroup of GL{d, M), { 1} C S" C 
GL{d,M). It defines an equivalence class of directional frames that we can call a partially 
framed fiag^^ The elements of our stabiliser S are given at the infinitesimal level by the 
formula (j62p . Therefore, an element s of can be written as 

g ^ gaiPo+a2P?+Q37'3+...+ad-l-Po"' (63") 

where ai, 02, 03, a^-i are constants. The element s of 5 acts on the coframe as, 

uj'' = s'jJ . (64) 

It is then easily checked that S is a commutative group of dimension d — 1. More precisely, 
in view of the fact that the various powers of the matrix Vo commute among themselves, 
one finds that the group composition of two elements of S is simply given by 

s{ai,a2, . ■ ■ ,ad-i) ° s(a'i,a2) • • • = + 0.1, <^2 + 0.2, ■ ■ ■ ,ad-i + a'd-i) ■ (65) 

[In mathematical terminology, S is a unipotent abelian subgroup of the Borel subgroup of 
GL{d,M).] Explicitly, the matrix elements of s (j63|) read as follow, 

a™ 

E ^Vn-l...V+^- (66) 

je(l,...,d-l)\3me{N*\mj=n-i} 

For instance in d = 3, the matrix s reads, 

/ 1 \ 

s= aiT^to^i 1 . (67) 

\ ilaj + a2)V^o)SV^o)\ ai7^(o)^2 1/ 
Therefore, in d = 3, the class (defined by the relation (I64p ) of coframes equivalent to some 
given coframe a;'* is explicitly given by: 



/2 2 I -n 2 1 



= uj"^ + alV^^,)'^2^'^ + {\oil + a2)V<^,)'^2V^a^'^iuj^ . (68) 



^'^Let us recall that a (complete) flag can be seen as the equivalence class of 'directional frames' with the 
following equivalence relations. A 'directional frame' given by the directions {v\, Vd} is equivalent to the 
'directional frames' {v'l, u^} constructed by picking a first direction along the vector «i (v'l oc vi), then a 
second direction v'2 belonging to 2-plane spanned by vi and V2 {v'2 (x V2 + avi), a third direction belonging 
to the 3-plane spanned by vi, V2 and U3 {v'j oc 113 + I3v2 + yvi), and so on up to a last direction along a 
vector v'„ which is an arbitrary vector in R''. In other words, the stabiliser of a flag is the full subgroup of 
lower triangular matrices of GL(d, R). 
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By duality between frames and coframes (< uj^,ej >= 5j), one can then easily deduce the 
corresponding equivalence classes of frames. The equivalence class of a frame ei is given by 



[ei] = {e[ = ejS-^^^\seS}. (69) 



For d = 3, we have explicitly, 

e[ = ei - aiP(o)^ie2 + {^aj - 02)^(0)^2^(0)^163 

e'2 = 62 - 01^(0)^263 

e's = e,. (70) 

Summary: The equivalence class of frames and coframes with respect to which the limiting 
values of M(o) and 'P(o) take the canonical values M^o) = ^ and 'P(o) = Vo (|60l) is described 
by the relations (j69p and (|64p which contain d — 1 arbitrary parameters. This equivalence 
class defines a privileged geometrical structure at the singularity, which we call a partially 
framed flag. For instance, in d = 3 this partially framed flag comprises: (i) one privileged 
direction 63, which is independent of the basic frame tj*, (ii) then the equivalence class of 
the vector 62 must lie in a privileged 2-plane containing 63, (iii) finally, after having chosen 
a representative 62 in the class [62] {i-e. ai is fixed), the third vector ei must lie in a privi- 
leged 2-plane defined by ei — Q!i7^(o)^i62 + |ai'P(o)^2'P{o)^i63 and 63. [This is different from 
a flag where ei would have been an arbitrary direction in the full R^.] 



It is remarkable that at the 'chaotic' BKL limit, we have a such geometrical structure left. 
At the singularity, we could have guessed that nothing is left from the metric structure 
because of the chaotic character of the asymptotic dynamics. Let us note that our results 
have been partially anticipated in [2] where a law of 'rotation of Kasner axes' was derived. 
This law is somewhat similar to our results ()68p but, actually, it has a different physical 
meaning. Indeed, the first Kasner axis I which is preserved in the law approximatively 
derived in ^ is supposed to belong to the 'growing' eigen axis Pi = Pi < 0, so that it would 
be a different axis that would be preserved during further collisions. One would then have 
no privileged direction at the singularity. As we have shown here, there exists a well defined 
geometrical structure at the singularity: a partially framed flag which is rather 'rigid' in the 
sense that it depends only on d — 1 arbitrary parameters (while a generic fiag would involve 
d{d — l)/2 arbitrary parameters). 



7 Conclusion 

In this paper, we started by reconsidering the asymptotic dynamics, in the vicinity of a 
spacelike singularity, of the fields for 'non-chaotic' Einstein-matter systems. We have out- 
lined a new proof that gives this asymptotic dynamics (which is essentially given, at each 
spatial point, by a monotone power-law solution in terms of the proper time). Our method 
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is based on the Iwasawa decomposition of the spatial metric and on the Hamiltonian for- 
mulation of Einstein-matter systems. As in references [141 I13j . we have used the Fuchs 
theorem to conclude. More precisely, we have defined an asymptotic system of equations 
which is a system of ODEs and we have also defined asymptotic constraints. Next, we have 
shown that the 'differenced variables' (i.e. the differences between the solution of the exact 
Einstein-matter constrained equations and the solution of the asymptotic Einstein-matter 
constrained equations) obey a Fuchsian system. A solution of the constrained asymptotic 
system, together with initial data, can thus be used to parametrize an exact constrained 
solution. The advantages of our formulation is that it is shorter, more transparent (the 
neglected terms are walls / subdominant walls) and that we avoid the problem of the sym- 
metry of the metric encountered in [lll|T3]. In appendix El we discuss the spatial domain 
on which a Fuchsian analysis can be applied in this context and point out that our method 
cannot be used in some zero-measure co-dimension 2 submanifolds (in reference |14[ I13j . 
there is a quite involved construction to deal with these submanifolds). We also showed 
that this problem originates in the 'spinorial' nature of the eigenvectors of the second fun- 
damental form around the submanifold where 2 (or more) eigenvalues coincide. 

We next turned to the 'chaotic' Einstein-matter systems and formulated a precise statement 
for the chaotic BKL behavior. This is achieved along the same lines as our formulation of 
the non-chaotic case. We parametrize, at each spatial point, the generic solution of the 
asymptotic behavior of the fields close to a spacelike singularity in terms of a constrained 
system of ODEs (and some initial data). Then we argue that the difference between the so- 
lution of the exact 'chaotic' constrained Einstein-matter system and the asymptotic system 
just defined satisfies a 'generalized Fuchs system'. We leave to others the task of proving 
that such 'generalized Fuchsian systems' admit a unique, asymptotical vanishing, solution. 
Our purpose here was mainly to formulate, in precise mathematical and physical terms, 
this asymptotic characterization of 'chaotic' solutions of Einstein-matter systems. 

Finally, we addressed the question of the existence of some asymptotic geometrical structure 
defined at the singularity for a chaotic system. We knew that some of the metric variables 
had finite limits at the singularity and it was therefore natural to wonder whether we could 
extract some geometrical structure from these limiting values. A first slight, we could 
expect that the chaotic nature of the asymptotic dynamics would destroy any structure at 
the singularity. We showed that it is not the case: partially framed flags can be defined at 
the singularity. These partially framed flags are, as their name indicate, more 'rigid' than 
flags and less 'rigid' than frames. 
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A Evolution equations of the Hamiltonian and momentum 
constraints 

The purpose of this appendix is to obtain the evolution equations for the exact Hamiltonian 
constraint Ti. (jlOp and exact momentum constraints (jlip in our gauge choices. Here we do 
not consider matter for simplicity and we work in a coordinate basis. These evolution equa- 
tions can be derived from the Bianchi identities and the evolution equations for the spatial 
metric. In this perspective, we have to know what are the evolution equations in the gauge 
we are interested in. 

The first order action of pure gravity m D = d + 1 spacetime dimensions vr*-^, N, N^] 

reads, 

S[gij,iT'^,N, N'] = j dx^d'^xioijT:'^ - NH - N'Hi) . (71) 

Its variation with respect to tt^^ can be understood as the definition of the vTjj in terms of 
gfjj, the variations with respect to gij, N, N'^ giving respectively the equations of motion, 
the Hamiltonian constraint and the momentum constraints!^ Note that the dot means a 
derivation with respect to x^. Let us now use the Einstein-Hilbert action, J d^x ^ — (^^^r'^'i? 
in the so-called Palatini formalism [i.e. the Christoffel symbols Typ are considered to be 
independent of the metric elements g^y) to determine the link between the Hamiltonian 
equations of motion 5S/5gij = and the usual Einstein equations '•^^G^i, = 0. The variation 

^■^We recall the following relationships, 
where fi = (0,i) and v — (0,j). 
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of the Einstein-Hilbert-Palatini action gives, 



where ^^^g' denotes the determinant of the spactime metric, '^'G^jy is the Einstein tensor and 
the second equaUty is obtained by integration by parts (we neglect the boundary terms). 
For simphcity, let us assume that A^* = on-shell (but keeping dN"^ ^ 0), and that most of 
the usual relations between T^^yp and the derivatives of g^^y are constrained to hold, namely 

T.0 -pO 1 9,1 , 

J- Oi — -L tO — o r "tT 

2 5 N 
r'oo = lN^g''g,,+g''NNjg 

^^jk = 2^^''i9ij,k + 9ik,j - 9jk,i) ■ (73) 

Moreover we require that T^jk is related to the T-^fco = T-'ofc via T^q^ = gg^'^T^ik- We can 
then verify that the terms in 5T can be ignored to compute the functional derivative of S 
with respect to gij, N and A^* (the coefficient in front of these variation in ()72p vanish). It 
is then straightforward to derive the following relations, 

^^lEL = Ngi-^'^^G'^ + N^g^^^G'^'^g'^) + 0{N'') (74) 

"9ij 

= 2<^'G°°iVV + 0(iV') (75) 
^ = -2iVrG^° + 0(iV'^). (76) 

Note that we did not write explicitly the terms proportional to Ni because we will work in 
the gauge Ni = 0. On the other hand, we have, 

-k'^ - + 0{N'') (77) 

dgij 

n (78) 



5S 



When identifying SSehp with 5S, one obtains 



6N 

-- -7f (79) 



n = -2g^N^ '">G°° + 0(iV'=) = - J-Goo + 0{N^) , (80) 
re = 2gN + 0{N^) = -^g'^Goj + 0{N^) . (81) 
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On the other hand, the equations of motion 5S/5gij = are found to be equivalent to the 
equations 

WQij ^ ]_'^gij J_ Q^^k^ _ (82) 

Note that our result ([82]) is linked to our choice of 'rescaled lapse' N as basic lapse variable. 
The result would be different if we were using the usual lapse A^. 

Let us now use the Bianchi identities := Vj^^^^C^^ = to derive the evolution equations 
for the constraints. We use the equality, 

^' ~ y3(DT^ 2^^"^ ^ ' 

When inserting in the expression of Bq the relationship between 7i, 7ii and the components 
of the Einstein tensor (jSO] [8T|) as well as the equations of motion ()82p . we obtain the 
evolution equation for 7^, 

1 N ■ ~ ■ 

—drH - -V'Hi - {ViN)W + 0{N^) = . (84) 
Ig z 

Note that 7ii = QijW and that it is a tensorial density of weight 1 while is a scalar density 
of weight -1. The covariant derivatives V^Jii and VjA^ must take into account these weights: 
for instance V^Tii = ViW = di{g^^7ij) (in a coordinate frame) and VjA^ = di{N^)/^. 
From i?j = 0, we get 

^(^-V.W) + O(iV^) = 0. (85) 

Note that TL \s a, scalar density of weight 2 so that ViH = gdiiTL/g). In our gauge choices, 
i.e. N = 1 and iV* = 0, the equations dMI) and ([SS]) read. 



drU = gV'ni + g,in\ (86) 
d^n^-ViH = 0. (87) 



B Fuchsian Systems 
B.l Fuchs Theorem 

The general form of a Fuchsian system ( \31\ \T7\ I35j and references therein) for a vector- 
valued unknown function u{t, x) = {ui{t, x), Uk{t, x)), defined on an open subset of RxM" 
with values in M^, is 

tdtu + A{x)u = t^'f{t,x,u,d^u), (88) 

where dxU denotes a finite number of derivatives of u with respect to the variables x (they 
are not restricted to be of first order); the function / is defined on (0, Tq] x Ui x U2 (where 
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Ui is an open subset of and U2 is an open subset of and takes values in W^; and 

fi > 0. ^ is an analytic k x k matrix- valued function defined on C/i. The system (I88p is 
said to be Fuchsian if the matrix A{x) and the function / fulfill the following conditions, 

• condition on A : the matrix A{x) is required to satisfy some lower boundedness 
condition. One sufficient condition, that has been used in several works |141 113j. is 
that there is a constant e such that Real(A) > e > 0, for each eigenvalue A of ^ at 
any point. Recently, this condition has been relaxed to requiring that there exists 
an < a < ^ such that |(t|-^o"° be bounded for a varying in the interval < o" < 1 
[35j . Essentially, this condition means that the real part of the eigenvalues of A are 
everywhere strictly larger than — /i. This is of particular interest for us, since the 
matrix A relevant in our case is nilpotent so that |ct|'^ grows like a power of logo" as 
<T — > (so that we can simply use any a in the interval < a < /i). 

• condition on f : the 'source term' / (after having factored t^) is required to be 
'regular', i.e. f must possess an analytic continuation in x, u and dxU and, as a 
function of t, must be continuous on [0,to] for some finite time to- 

For more precise conditions on A and /, we refer to |14j and references therein. 

Fuchs Theorem: If the system (j88|) satisfies the above conditions to be Fuchsian, then it 
possesses a unique solution u that vanishes as t — > +0 (see e.g [HIES])- 

Note : After the change of variable t = e~'^ (such that the singularity is now located at 
r +00), a Fuchsian system reads, 

drU — A{x)u = e~^'^ f{T,x,u,dxu), (89) 

where, essentially, / must be analytic in x, u, dxU and bounded in r as r ^ +00. This is 
the form we shall use in the text. 



B.2 Shift of the eigenvalues of A 

One may wonder if it is possible to have a more precise description of how fast the solution 
u of the Fuchsian system (j88p goes to zero when t goes to zero. To answer this question, 
let us rewrite the system ([88]) in terms of the variable u defined as follows, 

u = t~\ , < A < |U. (90) 

Inserting ([UU|) in the system ([55]) . it is straightforward to obtain 

tdtu + A{x)u = t^f{t,x,u,dxu) , (91) 
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where 

A = ^ + A1, /x = /x-A 

/ = f{t,x,t\t^d^u). 

If / is regular, then / is regular. Note that the eigenvalues of A have been shifted by 
A > compared to those of A, so that if A satisfies the lower boundedness conditions 
eigenvalues (^)> —ji, so does A (with the corresponding jl = fi — \). Therefore the 'shifted' 
system (j9ip is again Fuchsian and we know that it admits a unique solution u that vanishes 
when t — > +0. This tells us that the unique solution of the Fuchsian system ()88p that 
vanishes as i — > +0 actually vanishes as u = t'^u with u = o(l), i.e. as o{t^) for any A < /x. 

Summary : the 'shift' of u allows us to gain a more precise information about how the 
unique (asymptotically vanishing) solution u of the system ([88]) decays as t — > +0. If the 
r.h.s. of the Fuchsian system (j88p decays as t^/ then u is a 0{t^^'') for any e > 0. When 
using the r variable, this essentially means that a source term decaying as e"^"^ corresponds 
to a unique, asymptotically vanishing, solution decaying as e~^~'^ for any < /i_ < 

Note : For completeness, let us define the notations O and o. A function F{t,x,p) defined 
on (0, To] xUi X U2, where Ui, U2 are open subsets of M" and respectively, is said to be 
0{G{t)) if there is a constant C such that 

\F{t,x,p)\ < C\G{t)\ for t e (0,to], ix,p) G K. 

The notation F = o{G{t)) is used to indicate that F/G tends to zero uniformly on compact 
subsets of Ui X U2 as t ^ 0. 

C Subtleties occurring when some of the eigenvalues of the 
second fundamental form coincide in the non— chaotic case 

The usual AVTD approach uses a rather complicated construction to deal with the neigh- 
borhoods of points where some eigenvalues of the second fundamental form kij coincide 
[lll[13]. Such a complication is needed because the frame vectors that diagonalise kij with 
respect to gij are not analytic in x near such points. Here we consider the behavior of 
Iwasawa variables in these regions. For a full comparison one should carefully analyze the 
different slicing hypersurfaces in the two approaches: gaussian slicing = 1 in AVTD vs 
pseudo-gaussian slicing A^ = 1 in our case. 

As a simple example, let us consider gravity in D = 4 coupled to a dilaton. Let us consider 
for simplicity the (generic) case where two of the eigenvalues coincide on some submanifold. 
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We choose as one of the frame vectors, the (analytic) eigenvector 63 corresponding to the 
third eigenvalue (which is supposed to stay away from the other two). The two other 
analytic frame vectors are chosen to be orthogonal to 63 and to each other (they are linear 
combinations of the eigenvectors corresponding to the nearly degenerate eigenvalue). In 
this orthonormal basis (or dreibein), the coefficients of the metric are gab = Sab, while the 
coefficient of the second fundamental form are given by a matrix K which is of the following 
form, 



K 



/ a + c 
b 

\ 



b \ 

-a + c 
d J 



where a, b, c and d depend analytically on the spatial coordinates. The eigenvalues of K are 
c{x) lb -y/ a(x)^ + 6(x)^ and d{x). Therefore, two eigenvalues will coincide when 0(2;)^ + 6(x)^ 
vanishes, which means that both a{x) and b{x) must vanish. This happens generically on a 
line in the three dimensional space since it gives us two conditions a(x*) = and 6(x*) = 0. 
If we were in d spatial dimensions, the submanifold C where two eigenvalues coincide would 
again be defined by the vanishing of some a(x)^ + 6(x)^ and therefore be a codimension 2 
submanifold. For convenience, let us replace the quantities a(x*) and 6(x*) by p{x'^) and 0(x*) 
such that a(x*) = p(x*) cos 6'(x*) and 6(x*) = p{x^) sm6{x^). Since a and b have generically 
their values between -00 and +00, we have p G [0,oo] and 9 G [0,27r[. Let us consider the 



Kasner metric, expressed in the time r = — Int, g[o]{T,x^ 



-Kt 



[see [HKIS]). We have 



5[o] 



/ e'^'^ (cosh pT + cos 9 sinh pr) e^'^ sin 9 sinh pr 

e'^'^ sin 9 sinh pr e'^^ (cosh pr — cos 9 sinh pr) 

\ 








(92) 



To compute the Iwasawa variables corresponding to the metric (j92p . we use the following 
explicit formulas ([6]), 



3 



AA2 



1 1 , ,1 + cos ( 

- — cr In ( 

2 2^2 

1 1 , , 1 + cos ( 

- — cr H — In ( 

2 2^ 2 



1 — cos ( 
2 

1 — cos ( 



dT 



smt/le' 



pr 



-pr 



^1 + cos 9)eP'^ + (1 - cos 9)e- 



-pr 








-pr 



-pT\ 



(93) 



The crucial point is that the co-dimension 1 submanifold S defined by the equation 1 -|- 
cos 9 = plays a singular role in the formulas (j93p . Indeed, the exponential growing term eP'^ 
always appears in the combination (l-|-cos 9)eP'^ . Therefore, in the open domain l-|-cos 9^0 
and p = +Vcl'^ + b'^ > 0, we have the generic Iwasawa behavior that the AA's have a finite 
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limit as r — ^ +00, namely 

M\a,= lim Af\{r) 



1 + COS U 

while the /3's have the following asymptotic behavior, 

1 1 1 I + COS6' 

P ~ --(c + p)r- -ln( ), 

2 1 . ^ . 1 , , 1 + COS P 



/3 ~-o(c-p)T + -ln(: 



2^ ' / 2 ' 2 

(5^ = dT. (94) 

Note also that, because > 0, we have the usual asymptotic ordering < 0^. How- 
ever, we sec that the asymptotic limit N'^2(o), which depends only on spatial variables, 
becomes singular on S (codimension 1 submanifold (with boundary) where cos 9 = — 1) 
which is, in our d = 3 half-membrane ending on the line a = 6 = 0, where the 

eigenvalues of K coincide. More precisely, A/^^2(o) tends to +00 as ^ — tt"*" and tends to 
— 00 as — 7r~. Correlatively the behavior of the /3's become singular on E. We have the 
asymptotic behavior (5"-{t,x^) ~ p"(x')t + /?" (x*) where, e.g. (ilix^) = -^ln(^^t^^) and 
/3o (x*) = 5 In ( i+^°s^ ) both become singular on S, while pi = — ^(c + p), = —^{c — p) 
[we recall that p{x^) = +^/a?ix^^Y+^{x^]■ Note also that, if one sits on S, one has the 
asymptotic behavior J\f^2{T,x^) = and /3^(r,.T*) = — i(c — p)r, /3^(r,x*) = —^{c + p)t 
where the signs of the p terms in and 0^ arc exchanged compared to the asymptotic 
behavior outside of S. In particular, on S wc have, asymptotically, /3^(t, x*) > /3^(r, x*), 
which contrasts with the generic result that asymptotically < This unusual behavior 
is the sign that the coefficient of the e~^(^^~'^^) symmetry wall vanishes on S. As said above, 
in our treatment we neglected this possibility on the account that it is non generic (as the 
coefficient in question is a square). We see now that this non generic behavior necessar- 
ily occurs on some codimension 1 submanifolds ending on the codimension 2 submanifolds 
where 2 eigenvalues of K coincide. 

However, let us emphasize that the location of the singular codimension 1 submanifold S 
is not geometrically fixed, but is somewhat arbitrary apart from the fact that it necessarily 
ends on the codimension 2 submanifold C where 2 eigenvalues coincide. Indeed, let us 
show that, by using a suitable, a; -dependent local 50(2, M) transformation, one can move 
S around £, in a manner similar to an ordinary-life flag moving around its pole. More 
precisely, let us perform the following (spatially dependent) rotation of the first two vectors 
ei and 62 of our orthonormal frame. 







(5) 






























/ 







(95) 



One now finds 



,1 i ^ cos 9 siu (2o ) - siu 9 cos (2n ) ^ sin (2a - 9) 

' 1 + sin (2a) sin ^ + cos (2a) cos ^ 1 + cos (2a - ^) ' ^ ^ 
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which shows that the new singular surface S' corresponding to the rotated basis is now 
located at 6(x^) — 2a(x') = +tt. This shows that S is similar to a 'Dirac string' singularity: 
it is a gauge dependent singular submanifold, whose location can be shuffled around by 
using a gauge transformation. The non-analytic (actually singular) M — > oo, ... behavior 
of the Iwasawa variables on S is evidently problematic within a Fuchsian system approach, 
because it obliges us to work in an open region U of which does not contain S. The fact 
we just showed that the location of S can be moved around means that we can essentially 
bypass this technical problem by using simultaneously two separate Fuchsian systems, cor- 
responding to different choices of underlying frames (w*, Cj) in the analytic spatial manifold 
Mci, yielding finite values Iwasawa variables in two complementary open regions U and U'. 
Such a construction bypasses the analyticity problems near S and T,'. However, the codi- 
mension 2 submanifold C (where two eigenvalues coincide) remains excluded from these 
two complementary Fuch analyses. In other words, the present Iwasawa variables-based 
approach discussed here cannot cover in an analytic manner the measure-zero submani- 
folds of Md where 2 (or more generally n > 2) eigenvalues coincide. We leave this technical 
problem to further analyses. 



In this respect let us remark that the root of the problems linked, either in the AVTD 
or the Iwasawa approaches, to coinciding eigenvalues of kij admits a simple geometrical 
interpretation. The crucial point is that the eigenvectors of K, considered as functions 
of the auxiliary angle 6 introduced above (a = pcosO, b = psinO) depend on 6 in the 
same manner as a spinor would transform under a S0{3) rotation of angle 9. Indeed, the 
diagonalisation of the matrix K is easily checked to yield the following eigenvectors (with 
respect to the orthonormal basis ei, 62, 63 in which K is expressed). 



cos ~^ sin -62 with eigenvalue c + p = c + \/a'^ + b'^ , 



sin -ei — cos -62 with eigenvalue c — p = c — y/a?-\-l^ , 



V2 

V3 = 63 with eigenvalue d . (97) 

Note the appearance of the half angle 0/2 in vi and V2- This appearance means that if we 
follow the evolution of a diagonalising frame vi, V2, V3 along a closed loop in space around C 
(assuming the Jacobian of a{x) and b{x) never vanishes), the eigenvectors vi{x^) and V2{x^) 
will, upon their return to the same spatial point x (which corresponds to the same values 
of a(j;*) and 6(x*), but to an angle 6 = 2-k instead of 6* = 0), take final values opposite to 
their initial ones: 

vi{e = 2tt) = -vi{e = 0), V2{0 = 2tt) = -V2{9 = 0), Vi{e = 2tt) = Vi{9 = 0) . (98) 

This phenomenon explains why, in the Iwasawa approach, the presence of a 'line' C makes 
itself felt far away from C {i.e. on the 'singular' half-membrane S): indeed there is a 
non-trivial holonomy of Kasner frames around 



Evidently, one could discuss the richer case where there are several independent codimension 2 manifolds 
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D An illustrative Fuchsian toy model 



The aim of this appendix is to see, on a concrete example, the relationship between the 
structure of the source term and that of a solution of a Fuchsian system. The aim is not to 
sketch mathematical proofs (for which we refer to [14:\ [T3l [35] ) but to build some physical 
intuition. In this perspective, we study a toy model which is the most drastic simplification 
of the system (f55|) one could consider. We want to show that the unique solution of the 
system that goes to zero goes to zero like the source and not less fast than the source. To 
handle this, we use an iterative method. Finally, we also investigate a sligtly more involved 
toy model that is supposed to mimick the effect of the spatial gradients appearing in the 
system 



D.l A first toy model 

The first toy model we consider is the following, 

$-TT = 0, 

ir = e-^'"'^. (99) 

If we consider the asymptotic system to be 

/3[o] — '^[0] = , 

^[0] = 0. (100) 

The order zero solution is vrpj = v and /3[o] = vt + Po- Rewriting the equations ([99]) in terms 

of /3 = /3 — /3[o] and n = n — vrfoj gives 

/§-7f = 0, 

The idea is that (3 and vf go to zero as r ^ +cxd. We can try to solve the system by iteration 
assuming that f3 is small. The first iteration is obtained by replacing (3, in the source term 
of equations (|10ip . by its first order estimate which is zero, i.e. replacing e~'^^^^^ by 1: 

Ai] ~ ^[11 = ' 

The solution of this system is 7f[i] = ^sb^e^^ioCi-r+zSo) g^^^ ^ _J_^e-2«'(t'r+/3o) (^^^^^ 
that we did not add integration constants since we search solutions that go to zero when 
r — > +oo). The next step consists in replacing (3 in the source term by its first order 
estimate, 

P[2] ~ ^[2] = , 

^[2, = e-2«'(^^+^°)e-2"'(%l) ~ e-2"'(''^+^°)(l - 2u'(Ai])) • (103) 



of the type of £, together with more exceptional submanifolds where more eigenvalues coincide. 
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The solution of this system is given by vfp, = ^^e'^-C^+Z^^) - ^-^-^^s^^-^g-^^'^^m) and 
A,. = _L e-2»(-+'3o) _ ^2^^e-4-(/3[o]). Etc. 

Therefore, the solution that vanishes goes to zero like the source. We do not lose a e"^"^ and 
we conclude that the Fuchs theorem is too strong in this situation. For this very simple toy 
model, the exact solution of the system (I99p can be written as follows, 

*' (104) 



'2E - ie-2"'/5 

w 

where is a constant of integration and this integral gives explicitly, 



3 1 / e-2»"/3 

+ —^ln{l + Jl-^-—) + C, (105) 



V2E wV2E V 2Ew 

where C is a constant of integration. 

D.2 The effect of spatial gradients 

Let us now consider a less primitive toy model, mimicking walls with coefficients depending 
explicitly on the spatial derivatives of (3, 

$ = TT, 

n = a^(e-2'"^). (106) 

As asymptotic system, we take the same as in the previous example (jlOOp . The order zero 
solution is vrfoj = v{x) and /3[o] = v{x)t + f3o{x). Rewriting equations ()106p in terms of 
P = P — /3[o] and vr = vf — 7r[o] gives 

/5 = 7f, 

# = -2u;(9^OT + a^./3o + 9^^)e-2"'(^^+/5°)e-2'"(^) . 

The first iteration is again obtained by putting (3 is the left hand side to zero, 

vf[,j = e-2-(/5[oi)(%^ + ^) + 3fLre-^-^if^loi) (i07) 
w{v) w{v) 

is also of the form e~^'"^'^[''l^(a + 6r). Higher order iterations will give higher powers of 
r which are multiplied by increasing powers of the 'walls' e~^^[°l , see |36j for the structure 
of an all-order iterative example of such a Fuchsian system. Therefore, we can suspect that 
in this case, the solution decreases less quickly than the source by a polynomial in r. 
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